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GYROSCOPIC  CONTROL  AND  STABILIZATION 


Li-Sheng  Wang  *  P.S.  Krishnaprasad  f 

ABSTRACT.  In  this  paper,  we  consider  the  geometry  of  gyroscopic  systems  with 
symmetry,  starting  from  an  intrinsic  Lagrangian  viewpoint.  We  note  that  natural 
mechanical  systems  with  exogenous  forces  can  be  transformed  into  gyroscopic 
systems,  when  the  forces  are  determined  by  a  suitable  class  of  feedback  laws.  To 
assess  the  stability  of  relative  equilibria  in  the  resultant  feedback  systems,  we  extend 
the  energy-momentum  block-diagonalization  theorem  of  Simo,  Lewis,  Posbergh, 
and  Marsden  to  gyroscopic  systems  with  symmetry.  We  illustrate  the  main  ideas 
by  a  key  example  of  two  coupled  rigid  bodies  with  internal  rotors.  The  energy- 
momentum  method  yields  computationally  tractable  stability  criteria  in  this  and 
other  examples. 


1.  Introduction 

Geometric  control  theory  has  led  to  the  development  of  a  large  body  of  results  to  analyze 
and  design  nonlinear  feedback  systems.  There  is  a  beautiful  structure  theory  of  nonlinear  control 
systems  that  relates  internal  representations  (state  space  models)  to  external  representations  (input 
output  models).  Methods  from  the  differential  geometry  of  foliations  play  a  key  role  in  the  solution 
of  fundamental  problems  such  as  disturbance  decoupling,  non-interacting  control  etc.  cf.  [2i][33]. 
Inspired  in  part  by  the  success  of  the  linear  theory  and  partly  by  the  search  for  suitable  feedback 
invariants  [10],  the  paradigm  of  feedback  linearization  has  had  diverse  applications  [38].  Much  of 
this  work  is  however  concerned  with  generic  dynamics  and  often  does  not  specialize  well  to  the 
context  of  natural  mechanical  systems.  However,  there  is  a  growing  body  of  literature  devoted  to 
development  of  a  geometric  control  theory  for  hamiltonian  systems.  See  [12],  [33],  [22]  for  recent 
developments. 

This  work  was  supported  in  part  by  the  AFOSR  University  Research  Initiative  Program 
under  grant  AFOSR-  90-0105,  and  by  the  National  Science  Foundation’s  Engineering  Research 
Centers  Program:  NSFD  CDR  8803012. 

*  Institute  of  Applied  Mechanics,  National  Taiwan  University,  Taipei,  Taiwan,  R.O.C. 
j  Electrical  Engineering  Department  Sz  Systems  Research  Center,  University  of  Maryland. 

College  Park,  MD  20742. 
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A  parallel  intellectual  program  to  geometrize  mechanics  has  been  very  successful  and  has  had 
tremendous  impact  in  many  areas  of  mathematics  (e.g.  symplectic  geometry  and  topology)  and 
physics  (e.g.  gauge  theory,  geometric  quantization).  Some  of  the  roots  of  this  program  may  be 
traced  to  the  problems  of  stability  [4],  symmetry  and  reduction  [44]  [43]  [13]  [32],  and  investigation 
of  the  topology  of  phase  space  [42].  Modern  expositions  of  these  developments  may  be  found  in  [ij 
[5]  [6].  An  exciting  recent  result  is  the  block  diagonalization  theorem  for  simple  mechanical  systems 
with  symmetry  [39]  [41].  This  theorem  provides  refined  criteria  for  stability  assessment  for  natural 
mechanical  systems  by  careful  exploitation  of  the  underlying  geometric  structure. 

One  of  the  principal  goals  of  the  present  paper  is  to  demonstrate  the  applicability  of  geometric 
ideas  to  a  large  class  of  feedback  systems  derived  from  natural  mechanical  systems.  Here  the 
feedback  controls  are  of  the  gyroscopic  type.  A  key  point  is  that  simple  mechanical  systems  with 
symmetry,  when  subject  to  exogenous  forces  determined  by  suitable  classes  of  feedback  laws,  also 
admit  hamiltonian  and  lagrangian  structures.  In  [9]  [7]  the  hamiltonian  structures  so  derived  are 
viewed  as  deformations  (by  feedback  gains)  of  the  hamiltonian  structure  governing  the  open  loop 
unforced  system.  There,  the  methods  of  geometric  mechanics  such  as  reduction,  reconstruction 
phases,  and  the  energy- Casimir  algorithm  for  stability  analysis,  are  brought  to  bear  on  a  key 
example  of  rigid  body  control  using  external  torques  (as  implemented  by  gas  jets)  and  internal 
torques  (via  reaction  wheels/rotors),  and  the  relationships  between  these  two  methods  of  control. 
In  the  present  paper,  taking  an  intrinsic  lagrangian  viewpoint,  we  develop  a  systematic  theory 
of  gyroscopic  feedback  systems  with  symmetry.  Key  examples  of  such  systems  include  dual-spin 
satellites,  rigid  body  satellites  with  magnetic  torquers  etc. 

The  principal  reason  for  taking  a  lagrangian  viewpoint  is  that  it  leads  very  naturally 
to  the  incorporation  of  exogenous  forces/controls.  Furthermore,  in  the  setting  of  constrained 
nonholonomic  systems,  the  Lagrange-D ’Alembert  principle  is  the  basic  principle  of  modeling.  (See 
however  related  remarks  about  vakonomic  mechanics  and  the  role  of  variational  principles  for 
constrained  systems  in  [6].)  Good  representations  of  higher  order  tangent  bundles  together  with 
the  intrinsic/invariant  formulation  of  lagrangian  mechanics  lead  to  effective  modeling  of  the  systems 
of  interest.  Here  we  give  variational  principles  for  relative  equilibria  and  their  stability.  One  of 
the  contributions  of  this  paper  is  the  extension  of  the  (energy-momentum)  block-diagonalization 
theorem  to  gyroscopic  systems  with  symmetry.  Thus,  we  are  able  to  establish  a  set  of  refined 
stability  criteria  for  a  wide  class  of  feedback  systems  by  exploiting  fully  the  underlying  geometric 
and  group-theoretic  structures. 

The  outline  of  the  paper  is  as  follows.  Section  2  lists  relevant  notations  in  geometric  mechanics 
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and  group  action.  In  Section  3,  we  give  a  brief  treatment  of  lagrangian  mechanics  in  invariant 
form.  We  formulate  the  Lagrange-D ’Alembert  principle  in  geometric  terms.  Our  exposition  follows 
Vershik  and  Faddeev  [45]  in  large  part,  except  for  certain  conventions.  We  then  display  a  useful 
representation  of  the  second  tangent  and  cotangent  bundles  associated  to  S0(3).  In  Section  4, 
we  axiomatize  the  notion  of  gyroscopic  feedback  system  with  symmetry.  The  dual-spin  equations 
are  derived  as  an  example.  In  Section  5,  we  characterize  relative  equilibria  by  the  Principle  of 
Symmetric  Criticality.  We  give  a  careful  exposition  of  the  concept  of  relative  stability  in  the 
abstract  setting.  In  Section  6,  we  prove  the  block  diagonalization  theorem  for  gyroscopic  systems 
with  symmetry.  This  extends  naturally  the  previous  work  of  Simo,  Lewis,  Posbergh,  and  Marsden 
on  block-diagonalization  of  simple  mechanical  systems  with  symmetry.  A  distinctive  feature  of 
the  present  work  is  that  all  computations  are  done  on  the  velocity  phase  space  TQ ,  or  loosely, 
on  “lagrangian  side”.  We  are  aware  that  D.  Lewis  has  carried  out  a  similar  program  [27],  but 
the  present  work  was  done  independently  and  primarily  motivated  by  feedback  stability  problems*. 
The  correct  modification  of  the  amended  and  augmented  potentials  to  incorporate  gyroscopic  terms 
yields  stability  criteria  that  explicitly  display  said  terms.  This  is  further  made  clear  in  the  detailed 
example  of  two  coupled  rigid  bodies  with  internal  rotors  (the  multibody  dual-spin  problem  of  [48]) 
studied  in  Section  7. 

There  are  other  aspects  of  gyroscopic  feedback  controls  that  we  do  not  explore  in  this  paper 
but  we  think  are  quite  promising.  Control  strategies  based  on  bifurcation  of  relative  equilibria  may 
be  effective  in  a  variety  of  problems.  We  see  instances  in  [54],  and  [8].  In  the  present  setting  it  would 
be  worthwhile  to  investigate  bifurcations  with  respect  to  the  gyroscopic  feedback  parameters.  In  the 
context  of  the  dual-spin  problem,  this  has  been  carried  out  by  Ivrishnaprasad  and  Berenstein  who 
gave  a  bifurcation  diagram  [25].  Also,  in  [7],  the  authors  show  in  an  example  how  the  phenomenon 
of  geometric  phase  shift  is  affected  by  gyroscopic  parameters.  We  hope  to  discuss  these  aspects 
in  a  later  paper.  Some  of  the  results  in  this  paper  appeared  in  the  Pli.D.  dissertation  of  Li-Sheng 
Wang  [46]. 

We  would  like  to  thank  Professor  Jerrold  Marsden  for  critical  and  helpful  advice  on  an  earlier 
version  of  this  paper. 


*  A  feedback  law  is  a  rule  that  determines  exogenous  forces/controls  as  functions  of  the  current 
state  of  a  dynamical  system.  The  feedback  stabilization  problem  is  to  find  feedback  laws  that 
achieve  prescribed  characteristics  in  the  closed-loop  system. 
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2.  Notations  and  Preliminary 

In  this  section,  we  introduce  the  notations  used  in  later  sections.  First,  we  collect  together 
some  basic  geometric  objects  of  lagrangian  mechanics.  We  follow  Vershik  and  Faddeev  [45]  closely, 
filling  in  details  where  needed. 

Let  Q  be  a  smooth  manifold  with  local  coordinates  x,  TQ  be  the  tangent  bundle  of  Q  with 
local  coordinates  (x,v),  and  7 r  :  TQ  —>  Q  be  the  canonical  projection  from  TQ  to  Q  .  Let  TTQ  be 
the  second  tangent  bundle  with  local  coordinates  (x,  v,  u,  w).  Let  T(x,v)  be  the  tangent  space  of 
TQ  at  (.r,  u),i.e.  T(X<V)TQ  ,  and  denote  the  vertical  tangent  subspace  of  T(x,v)  consisting  of  vectors 
tangent  to  the  fiber  of  TQ  by  T^x  ^ .  In  local  coordinates,  each  vector  in  can  be  written  as 

(0,  w),  for  some  w  £  TXQ .  The  tangent  map  Tbr  of  the  canonical  projection,  TV  :  TTQ  — >  TQ , 
can  be  expressed  in  local  coordinates  as  Ttt^x^v)(u,  w)  —  u  6  TXQ,  which  projects  a  vector  in  T(XtV) 

to  its  horizontal  component.  Define  the  map  7(X)U)  :  TXQ  — >  T(X<V)TQ\  u  (0,  u)  £  dis¬ 

establishes  an  isomorphism  between  TXQ  and  2)^  ^ .  Let  X  :  TQ  —*  TTQ  be  a  vector  field  on 
TQ  .  X  is  called  a  vertical  vector  field  if  X(x ,  v)  £  Ty.  ^  or  X(x,  v)  =  (0,  w) ,  for  some  w  £  TXQ  . 
This  is  equivalent  to  saying  Tix  ■  X  —  0.  There  is  a  unique  vertical  vector  field  Xp^  called  the 
principal  vertical  field  defined  by  Xpv(x,  v )  =  ~i(x,v)  •  v  =  (0,  v). 

* 

Now  we  consider  the  dual  spaces.  Let  w  be  a  1-form  on  TQ .  It  is  said  to  be  horizontal  if 
for  all  vertical  vector  fields  X,  oj(X)  =  0,  or  in  local  coordinates,  u(x,  v )  =  (a,  0),  for  some 
a  £  TfQ.  The  dual  of  the  map  Tiv(XtV)  :  T{x,v)TQ  —  TXQ  can  be  defined  implicitly  through 

T*{x,v)  :  T;Q  -  T{^V)TQ, 

(TTC(x,v)a>  (u>w))  =  (°b  T~(x,v)  (u,w))  =  ( a,u ), 

where  a  £  TfQ .  Thus,  in  local  coordinates,  =  (cr,  0) .  Similarly,  the  dual  of  7(2jV), 

denoted  by  7^  :  T^x  v)TQ  ->  TfQ ,  is  defined  by,  in  local  coordinates, 

(l{x,v)(ai  P)i  u)  =  ((a>  P)i  7 (®.u)  '  u )  —  (Pi  u )• 


Equivalently,  7(X)V)(ai  (5)  =  (5  ■ 

With  these  dual  mappings,  we  define  the  bundle  map  r  :  T*{TQ) 
t[x,v)  =  T**X,V)  ‘  7(*r,«) •  In  particular,  for  (a,/?)  £  T*xv)TQ  ,  we  have 

T(x,v)  («,  P)  =  Tir(x,v)  ■  l{x,v)  (a>  P)  ~  (Pi  °)' 


T*(TQ)  as 


(2.1) 


Thus  T(X,V)  maps  any  cotangent  vector  ( covector )  to  a  horizontal  covector.  Globally  r  maps  any 
1-form  on  TQ  to  a  horizontal  1-form  on  TQ  .  On  the  other  hand,  we  may  define  a  bundle  map  from 
the  second  tangent  bundle  into  itself,  r»  :  TTQ  — >  TTQ  as  t^X}V)  —  '7(x,v)  ■  l°ca^ 
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coordinates  we  can  associate  to  each  !u,w)  £  T,~  , \T0  t  ,  Jv  i  T  ^ 

v  )  v?  r*(x,v)(^?  w)  —  (0.,  u).  In  other  words, 

maps  any  second  tangent  vector  to  a  vertical  tangent  vector,  and,  globally  r*  maps  a  vector 

field  on  TQ  to  a  vertical  vector  field  on  TQ  . 

To  treat  second-order  equations,  we  need  the  following  concept.  A  vector  field  on  TQ , 

X  G  * (TQ)  ’  iS  a  Spedal  vector  field  if  only  if  rmX  =  XpV .  In  local  coordinates,  assuming 
X(x,v)  =  (u,w),  it  says  t*X{x,v)  =  (0,  »),  which  is  equivalent  to  the  condition  u  =  v.  It 
then  follows  that  this  definition  of  special  vector  field  X  is  the  same  as  saying  X  gives  rise  to  a 
second-order  equation  on  Q,  cf.  p.  213  of  Abraham  and  Marsden  [i], 

Let  T(x,v)  denote  the  space  of  horizontal  covectors  at  (x,  v)  in  TQ .  Define  the  map 
:  T(Z)  to  be,  in  local  coordinates,  a(x>v){a,  0)  =  a,  for  a  £  TfQ.  This  map  will 

be  used  later  in  defining  the  Legendre  transformation.  The  maps  7,  r,  r,,  and  a  are  intrinsic  and 
do  not  depend  on  choices  of  local  trivializations  of  the  bundles  involved. 

Now  we  collect  together  basic  notions  of  group  actions  on  riemannian  manifolds  necessary 
to  discuss  natural  mechanical  systems  with  symmetry.  Let  (Q,<  »  be  a  manifold  with 

riemannian  metric  <•,•>.  We  sometimes  write  K(x){vx,  wx)  =«  vx,  wx  >x,  for  xeQ,  and 

vx,Wx  G  TXQ .  The  liemannian  metric  induces  a  vector  bundle  isomorphism  K **  :  TQ  — *  T*Q 
defined  by 

(Ii\vx),  wx)x  =  <  vx,  wx  »*,  for  all  vx,  wx  £  TXQ, 


where  {•,  -)x  denotes  the  pairing  between  elements  in  TfQ  and  TXQ  .  Here,  and  in  what  follows, 
the  notation  (•,  •)  is  used  to  denote  the  dual  pairing  between  appropriate  spaces.  By  the  Biesz 
Representation  Theorem,  this  isomorphism  is  well  defined  and  we  may  write  10  =  (I0)~ 1  : 
T  Q  -+TQ,  which  is  also  a  fiber-preserving  mapping.  By  definition,  for  ax  £  T*Q  and  wx  £  TQ , 

we  have  (ax,  wx)x  —  <C  10  •  olx,  wx  ~^>x.  Via  the  isomorphism  10  ,  an  induced  inner  product  on 
T*Q  can  be  defined,  for  ax,  f3x  £  T*Q , 

<  OLx,  (3X  >T.Q  =  <  10  .  ax,  10  •  f3x  .  (2.2) 


Let  G  be  a  Lie  group,  and  §:G  X  Q  Q  be  a  group  action  of  G  on  the  manifold  Q  . 
We  shall  nse  the  notations  x)  =  $<,(»)  =  g  ■  x  interchangeably  to  denote  this  action.  The 
tangent  lift  associated  with  d>  is  defined  as  <5  J  =  TT,  :  TQ  ~  TQ ,  or,  in  local  coordinates, 
$5(.r,  v)  —  ($g(x),  Tx<t>g  ■  v).  The  cotangent  lift  <I>r  associated  to  i>  on  the  cotangent  bundle 
T*Q  ,  <hT  :  G  X  T*Q  ->  T*Q ,  is  (a*)  =  T*$g-x  •  ar,  where  1  is  the  dual  of  T$g-i .  In 

local  coordinates,  we  have  (<J>f  (*,  a),  (g  -  ar,  v))g.x  =  <a,  Tg.x<hg~i  •  v)x.  It  is  straightforward  to 
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verify  that  the  tangent  lift  and  cotangent  lift  are  both  well-defined  actions  on  the  spaces  TQ  and 
T*Q  respectively. 

Let  the  Lie  algebra  of  a  Lie  group  G  be  denoted  by  Q ,  with  its  dual  Q* .  Recall  that  the  Lie 
algebra  Q  is  identified  as  the  tangent  space  to  G  at  the  identity  element  e  or,  equivalently,  the  set 
of  left  invariant  vector  fields  on  G ,  cf.  also  [34].  Given  £  £  Q ,  for  a  group  action  $  on  a  manifold 
Q ,  we  define 

=  ^expf^)  E  TXQ, 

a€  6-0 

the  infinitesimal  generator  of  the  action  corresponding  to  £ .  The  group  G  acts  on  Q  through  the 
adjoint  action 

Ad  :  G  X  Q  -  0;  (g,  £)  ^  Te(Rg- 1  0  Xs)£  =  Adff£,  (2.3) 

where  Lg ,  Rg  denote  the  left  and  right  translation  of  a  group  element  by  g  £  G ,  respectively.  The 
map  g  Adg  is  also  called  the  adjoint  representation  of  G  in  Q .  The  infinitesimal  generator  of 
this  adjoint  action, 

£C'(r?)  =  XT  Adexpe^(rj) 
ae  e  =0 

can  be  shown  to  be  equal  to  the  Lie  bracket  of  £  and  g ,  namely, 

Zoiri)  =  [  £,  V  ]  =  ad^g. 

(We  follow  the  sign  convention  for  Lie  brackets  used  in  [1].)  The  group  G  also  acts  on  the  dual  of 
the  Lie  algebra  G*  through  the  coadjoint  action 

Ad*  :  G  x  Q*  -  £T;  (g,  p)  .*-»■  Adj-iju, 

which  is  defined  by,  ( £)  =  (//,  Ada£) ,  for  all  £  G  17 ,  where  (•;,  •}  denotes  the  duality  pairing 
on  Q*  x  Q .  The  corresponding  infinitesimal  generator,  £5*  can  be  shown  to  be  determined  by 

{£<?«(/£),  V)  =  -(id  [£,  ’?])  =  -  (adj/r,  t?), 

for  all  g  G  £7.  From  the  previous  definitions,  it  is  easy  to  establish  the  identity,  (ad£/i,  77)  = 
(/z,  ad^g).  With  the  above  structure,  group  G  acts  on  Q  and  Q*  through  the  actions  $  and  Ad* 
respectively.  A  map  J  :  Q  — *  Q*  is  called  Ad*  -equivariant  if  J  o  $<,  =  Ad*- 1  0  J.  For  p  E  G* ,  we 
define  the  isotropy  subgroup  of  p  by 

Gfj,  =  {  g  E  G  :  Ad*p  =  /t  },  (2.4) 

with  isotropy  Lie  algebra 
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Gy.  =  {  V  e  Q  ad*n  =  0  }, 


(2.5) 


which  can  be  shown  to  be  a  subalgebra  of  Q . 

The  notion  of  invariance  is  now  ready  to  be  introduced.  A  riemannian  metric  is  G-invariant 
if  it  is  invariant  under  the  pull-back  of  the  mapping  <f>3,  i.e.  for  all  g  6  G,  d>*  •  K  —  K ,  or 
in  local  coordinates,  K(x)(v,  w)  =  K(g  ■  x)(Tx$g  ■  v,  Tx$g  •  w),  Vv,  w  E  TXQ.  It  then  follows 
that  the  inner  product  on  T*Q  defined  in  (2.2)  is  invariant  under  the  cotangent  lift,  namely, 
<  o-x,  fix  >T'Q  —  <  axi  Px  >Tm Q  t  all  9  £  G .  This  can  be  shown  by  using  the  following 
identities, 

Ii“  •  =  Tx$g  ■  A3  •  ax,  for  a*  £  T*Q,  (2.6a) 

Kh  •  Tx$g  •  wx  =  •  Iib  •  wx,  for  wx  6  TQ.  (2.66) 

Similarly,  a  1-form  Y  G  T*Q  is  called  G-invariant  if  -Y  =  Y ,  for  all  g  G  G .  A  smooth  function 
V  :  Q  — IR  is  a  G-invariant  function  on  the  manifold  if,  for  all  g  G  G ,  V  ($g(x))  =  V(x).  A  vector 
field  Y  on  Q  is  a  G-invariant  vector  field  if  for  all  g  G  G ,  ($s)*-y  =  Y ,  or  Y(x)  =  T^g-Y(g~1  -x) , 
for  x  e  Q ,  g  e  G . 

Recall  that  a  differential  operator  on  the  full  tensor  algebra  can  be  defined  from  its  restrictions 
on  functions  and  vector  fields,  cf.  the  theorem  of  Willmore  [52].  Accordingly,  the  Lie  derivative  of 
a  vector  field  on  the  tensor  algebra  can  be  found  from  its  Lie  derivative  on  functions  (directional 
derivative)  and  Lie  derivative  on  vector  fields  (Lie  bracket).  The  following  two  lemmas  are  essential 
to  the  developments  in  Section  6.  Their  proofs  can  be  found  in,  e.g.  [1],  and  [46]. 

LEMMA  2.1 

Let  Y ,  K  be  G-invariant  vector  field  and  riemannian  metric,  respectively.  Their  Lie 
derivatives  with  respect  to  vector  field  tjq  ,  where  77  G  G,  vanish,  namely,  LVqY  =  0,  and 

lVq  k  =  0. 

I 

LEMMA  2.2 

Let  AT,  A 2  G  %(Q),  the  vector  fields  on  Q  .  For  g  £  Q ,  we  have 

Lvq  <  Al5  A2  »x  =  «  LVq  AT,  A2  »x  +  <  AT,  Lvq  A2  »*  • 
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3.  Lagrangian  Mechanics  in  Invariant  Form 

Lagrangian  mechanics  provides  a  systematic  formulation  of  mechanical  problems  from  a 
unified  point  of  view.  In  contrast  with  working  on  the  cotangent  bundle  as  in  most  of  hamiltonian 
mechanics,  lagrangian  mechanics  formulates  the  problems  on  the  tangent  bundle,  or  the  velocity 
phase  space.  As  we  shall  see,  it  admits  greater  freedom  in  interpreting  and  formulating  intuitive 
physical  notions  such  as  exogenous  forces  and  the  principle  of  virtual  power.  In  this  Section,  we 
present  the  invariant  form  of  lagrangian  mechanics  through  local  representations  and  show  that  the 
invariant  form  of  the  Lagrange-d’Alembert  Principle  gives  rise  to  the  Euler- Lagrange  equations  in 
local  coordinates.  Moreover,  a  similar  equation  can  be  applied  to  model  mechanical  systems  on  the 
special  orthogonal  group  (rotation  group)  50(3)  in  a  global  sense. 

Let  Q  be  a  smooth  manifold  viewed  as  the  configuration  space.  Let  TQ ,  TTQ  denote  its 
tangent  bundle,  second  tangent  bundle,  respectively.  Since  lagrangian  mechanics  is  about  second- 
order  equations,  we  need  to  consider  the  corresponding  elements  in  the  jet  spaces  of  Q  ,  namely  the 
second  tangent  vectors.  Let  L  :  TQ  — »  IR  be  a  smooth  function  (or  Lagrangian).  The  corresponding 
differential  1-form  dL  :  TQ  -+  T*{TQ)  can  be  written  in  local  coordinates  as, 

(dL(x,v),  (u,  w))  -  TL{XtV)  ■  (u,  w)  (3.1) 

or,  in  terms  of  Frechet  partial  derivatives,  dL(x,v)  =  (D\L(x,  v),  D2L(x,  v)) .  The  horizontal 
1-form  0 on  TQ  corresponding  to  L  is  defined  to  be,  cf.  (2.1), 

0L  ■=  r-dL.  (3.2 a) 

In  local  coordinates,  0L(a;,  v)  =  r(x>t,)  •  dL(x,v )  =  (D2L(x,v),  0) .  Taking  the  exterior  derivative 
of  the  1-form  0^ ,  we  associate  to  L  a  2-form  on  TQ ,  :  TTQ  x  TTQ  — >  IR,  defined  as 

nL  =  —  d  Ql.  (3.3 a) 

If  we  write  0l  =  DV,L  dx\  then,  by  taking  exterior  derivatives  on  both  sides,  we  get 

LIl  =  -dDv.L  A  dx{  =  DxjDv,L  dx{  A  dxj  +  DviDv,Ldxi  A  dvj .  (3.3 b) 

On  the  other  hand,  let  («i,tt?i),  (tt2,ti72)  6  T(x,v)TQ  .  From  (3.3a),  we  derive  the  following  formula 
in  local  coordinates, 

Q,l(x,v)((ui,  w i),  (uo,  w2))  =  (D\D2L(x,v)  ■  u2)  ■  u-i  +  (D2D2L(x,v)  •  w2)  •  u\ 

(3.3c) 

—  {D\D2L(x,  n)  •  tii)  •  ^2  —  {D2D2L{x,  v)  ■  wQ  ■  u2 
Next,  we  give  the  intrinsic  form  of  Legendre  transformation  which  maps  the  velocity  phase 
space  to  the  momentum  phase  space.  The  Legendre  transformation  corresponding  to  the  Lagrangian 
L  can  be  defined  as 
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lL  :  TQ  -*  T*Q;  (a:,  v)  h->  (a:,  cr(Xi„)  •  0L(x,  vj)  ,  (3.4) 

or,  equivalently,  Il(x,  v )  =  (a;,  D2L(x,v)) ,  cf.  the  definition  through  fiber  derivatives  in  pp.  209, 
219  of  [i]. 

Assuming  now  that  £l  is  a  diffeomorphism  (or  L  is  hyperregular),  we  have  l~^  :T*Q  -+  TQ. 
(This  condition  implies  that,  in  local  coordinates,  D2D2L(x,v)  is  nonsingular.)  Denote  the  space 
of  &-forms  on  a  manifold  M  as  zok(M).  By  the  pull-back  of  ij^- ,  (f!£  lY-.za\TQ)-*vo\T*Q), 
we  can  define  a  2-form  on  T*Q  as  w0  =  Although  f lj_,  is  L -dependent,  u0  defined 

above  is  invariant  under  the  change  of  L .  In  fact,  letting  ( x,p )  be  local  coordinates  of  T*Q  ,  where 
( p ,  w)  —  (D2L(x,v),  w) ,  for  w  E  TXQ  ,  it  can  be  shown  that  uq  —  dx  A  dp,  which  is  the  canonical 
svmplectic  2-form  on  the  cotangent  bundle. 

Thus,  when  the  Legendre  transformation  is  diffeomorphic,  the  two  approaches,  either  based 
on  the  cotangent  bundle  or  directly  on  the  tangent  bundle,  are  equivalent.  Moreover,  ( T*Q ,  u0)) 
(TQ,  f II)  are  both  symplectic  manifolds,  carrying  associated  Poisson  structures. 

REMARK  3.1 

The  closed  2-form  Ql  in  (3.3a)  is  well-defined  for  every  Lagrangian  L.  It  is  however 
nondegenerate,  and  therefore  a  symplectic  structure,  only  when  L  is  regular.  For  a  singular  or 
irregular  L,  fit  becomes  presymplectic ,  namely  is  no  longer  of  maximal  rank.  Discussions  of 

this  case  may  be  found  in,  e.g.  [14]  [15]. 

1 

With  the  symplectic  2-form.  one  constructs  a  correspondence  between  vector  fields  and 
1-forms,  IIl  :  wx(TQ)  — »  %(TQ )  through,  for  u  £  zv  1(TQ), 

fii(ni(w),  Z)  =  u  (  Z  ),  VZe%(TQ).  (3.5) 

In  terms  of  the  inverse  of  Ilr, ,  an  alternative  expression  is  Ql(X,  Z)  =  II£1(A)(Z),  for  all 
Z  G  %(TQ),ox 

n i\x)  (•)  =  nL  ( x,  •)•  (3.6) 

It  can  be  shown  that  II n  maps  horizontal  1-forms  to  vertical  vector  fields  [45]  [46].  Now  define  the 
energy  function  on  TQ  ,  Hl  :  TQ  — >  1R,  as, 

II L  =  dL  (Xpv)  -  L ,  (3.7) 

where  Xpv  is  the  principal  vertical  field  defined  in  Section  2.  In  local  coordinates,  we  have 
Hl(x,  v)  =  (£l(x,v),  v)  -  L( x,  v )  which  is  exactly  the  same  notion  as  the  energy  defined  in  p.213 
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iii  [1].  In  particular,  tlie  function  dL(Xp]/ )  is  sometimes  called  the  action  corresponding  to  L. 
From  the  energy  function  HL  on  the  velocity  phase  space,  define  the  hamiltonian  on  the  momentum 
phase  space  as 

H  :  T*Q  -  1R;  II  =  HL  o  TL\  (3.8) 

The  hamiltonian  system  (T*Q,uq,  H)  is  the  customary  object  of  study  in  hamiltonian  mechanics. 
The  Lagrangian  vector  field  determined  by  L  is  defined  as, 

XHl  =  nL(  dHL  ),  (3.9a) 

or,  equivalently,  f Il(Xhl,  Z )  =  dHi(Z),  for  all  Z  6  %(TQ).  In  local  coordinates,  the  matrix 
form  is  X Z  —  VHlT  Z.  Thus  we  may  write  the  Lagrangian  vector  field  as 

XHl  =  ([fiL]-1)TVifL.  (3.96) 

One  could  also  think  of  Xhl  as  the  Hamiltonian  vector  field  corresponding  to  the  Hamiltonian  Hl 
on  the  symplectic  manifold  (TQ,Hl),  and  thus  Hi  is  a  first  integral  (conserved  quantity)  along 
the  vector  field  XuL  .  We  say  that  we  can  define  consistent  equations  of  motion  if  such  an  X 'hl 
exists.  It  can  be  shown  that  Xhl  is  a  special  vector  field  and  thus  gives  rise  to  a  second-order 
equation. 

Now  we  introduce  the  important  notion  of  Lagrangian  force.  Recall  that  in  lagrangian 
mechanics  [45],  virtual  displacements  can  be  thought  as  special  vector  fields  on  TQ ,  and  forces 
can  be  modeled  as  horizontal  1-forms  on  TQ  .  For  a  Lagrangian  L  ,  the  associated  Lagrangian  force 
on  a  virtual  displacement  A',  Fl(X),  is  defined  through 

Fl(X)(Z)  =  flL(X,Z )  -  dHL{Z),  VZG  3e(TQ)  (3.10a) 

The  Lagrangian  force  Fl(X )  is  a  1-form  on  TQ.  This  1-form  can  be  shown  to  be  a  horizontal 
1-form  on  TQ  .  In  fact,  in  local  coordinates,  with  X(x,  v )  =  ( v ,  w ) ,  we  have 

Fl( X)(x,v)  (u,w2)  =  (— D-lD-2L(x,v)  •  v  —  DoDoL^x,  u)  •  w  +  D\L(x,v)  )  •  u.  (3.106) 

Thus  it  is  a  well-defined  force. 

Definition  in  (3.10a)  holds  even  for  L  singular,  cf.  Remark  3.1.  If  L  is  hyperregular,  we  may 
write,  cf.  (3.5),  Fl(X )  =  II^1(A)  -  IIIl.  This  is  the  definition  used  in  [45]  for  the  Lagrangian 
force.  In  the  above  setting,  the  Lagrange-d’Alembert  Principle  can  be  now  stated  in  the  following 
form. 
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LAGRANGE  D’ALEMBERT  PRINCIPLE  3.2 

For  a  holonomic  mechanical  system,  on  the  virtual  displacement  (special  vector  field)  that 
determines  the  real  trajectory  of  motion,  the  sum  of  the  Lagrangian  force  and  the  exterior  force 
vanishes. 

H 

For  natural  systems,  the  Lagrangian  force  consists  of  resultant  force  of  inertia  and  forces 
coming  from  the  potential  energy.  Thus  the  principle  here  corresponds  to  the  classical  d’Alembert 
principle,  see  e.g.  [26].  As  discussed  in  [26],  a  fundamental  entity  in  analytical  mechanics  is  virtual 
work,  instead  of  the  classical  notion  of  force.  Here  we  present  a  unified  treatment  in  terms  of 
horizontal  1-forms,  where  the  classical  forces  are  represented  by  the  coordinates  of  this  1-form. 

Let  u;  be  an  exterior  force  or  a  horizontal  1-form.  The  d’Alembert  principle  says  that 

Fl(X)  +  a;  =  0,  (3.11) 

where  is  a  special  vector  field.  The  trajectories  of  motion  of  the  mechanical  system  with 
Lagrangian  L  obey  the  flow  of  this  vector  field.  In  the  absence  of  any  exterior  force  and  L  being 
regular,  from  (3.10a),  we  write  9,L(X,Z)  =  dHL(Z),  for  all  Z  6  %(TQ),  which,  by  definition  of 
Xhl  ,  implies  that  X  =  X hl,  i.e.  the  Lagrangian  vector  field  gives  the  real  trajectories  of  motion.- 
Now  we  express  the  d’Alembert  principle  in  local  coordinates  where  (3.11)  reads 

Fl( X)  (x,  v )  +  u(x,  v )  =  0. 

Letting  u  =  (a,  0),  X(x,v)  =  (v,  w),  we  have,  cf.  (3.10b), 

~D1D2F(x,v)  ■  v  -  £>2F>2 l(x,v)  ■  w  +  DiL(x,v)Sj  ■  u  +  a  ■  u  =  0,  V  u  6  TXQ. 

By  including  time  derivatives  as  v  =  x ,  w  =  v ,  we  get 

^D2L(x,v)-u  =  D\L{x,v)-u  +  a  ■  u,  \/ueTxQ.  (3.12) 

Integrating  both  sides  with  respect  to  the  variable  t ,  this  equation  can  be  rewritten  as 

rv  /p  'J' 

DnL(x,v)  •  ut  dt  =  /  ( DiL(x,v)-u  +  a-  u)dt. 

Jo 

This  corresponds  to  the  Principle  of  Virtual  Power  in  analytical  mechanics,  cf.  e.g.  [51]-  The 
tangent  vector  u  is  sometimes  called  test  function.  In  the  case  that  the  pairing  is  nondegenerate, 
for  example,  in  the  finite  dimensional  case,  w-e  can  write  (3.12)  as 

jtD2L{x,v)  =  DiL(  x,v)  +  a ,  (3.13) 

which  is  the  classical  form  of  the  Euler-Lagrange  equation. 


D2L{x,v)-u  —  / 

0  Jo 
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EXAMPLE  3.3.  (On  the  group  50(3)) 

Now  we  illustrate  the  Lagrange-d’Alembert  Principle  in  the  setting  of  the  Special  Orthogonal 
Group  50(3)  as  configuration  space.  Recall  that  each  element  A  in  50(3)  is  an  element  in 
05(3),  the  group  of  all  3  X  3  nonsingular  matrices,  which  satisfies  the  condition  AT A  =  1  and 
det(A)  =  1 .  Let  the  operator  **  denote  the  natural  isomorphism  from  ]R3  to  so{ 3) ,  the  space  of 
3x3  skew-symmetric  matrices,  defined  by 

f  W\\  (  0  ~™3  W-2  \ 

I  W2  )  =  I  W3  0  -Wi  .  (3.14) 

\W3  J  \-Wo  Wi  0  J 

Given  A  £  50(3),  recall  that  (A,  AO)  is  an  element  in  550( 3).  In  rigid  body  mechanics,  the 
variable  ft  corresponds  to  the  instantaneous  angular  velocity  of  the  motion  in  body  coordinates. 
With  the  representation  of  the  elements  in  the  second  tangent  bundle  TTSO( 3)  as 

^A,  AO,  An,  A(nO  +  w)^j  ,  (3.15) 

and  the  trace  pairing  in  05(3) 

(A,  5)  =  Ur(ATB),  for  A,  5  £  05(3),  (3.16) 

we  can  write  elements  in  5*550( 3)  canonically  as 

(A,  AO,  A(60  +  a),  Abj  .  (3.17) 

Here  u,  w ,  a,  b  are  vectors  in  IR3 ,  and  the  pairing  between  5*550(3)  and  5550(3)  becomes 

(^A,  AO,  A(60  +  a),  Ab^j  ,  ^A,  AO,  Au,  A(nO  +  w)^ )  —  a-u  +  b-w. 

We  remark  here  that  these  parametrizations  of  5550(3)  and  5*550(3)  are  globally  defined 
via  the  embedding  of  50(3)  in  G5( 3).  Our  goal  has  been  to  make  the  pairing  analogous  to  that 
on  Euclidean  space.  The  global  representations  (3.15),  (3.17)  of  the  second  tangent  bundle  and  the 
dual  of  the  second  tangent  bundle  on  50(3)  also  prove  to  be  useful  in  computing  the  derivatives  or 
variations  of  a  function  (Lagrangian)  on  550(3)  and  in  deriving  the  reduced  Poisson  bracket  [49]. 
In  the  following,  we  state  the  Lagrange-d’Alembert  Principle  in  terms  of  these  representations. 

On  550(3),  let  a  system  be  described  by  a  Lagrangian  5.  The  Lagrange-d’Alembert 
Principle  in  the  invariant  form  (3.11)  applied  to  motions  on  50(3)  gives  rise  to  the  Euler-Lagrange 
equation,  namely  for  all  Au  £  5^50(3), 

525(A,  Aft),  Au)  =  (D\L(A,  Aft),  Au)  +  (a,  Au),  (3.18) 
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where  a  is  the  exterior  force.  See  [46]  for  detailed  discussions. 
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4.  Gyroscopic  Control 

In  this  section  we  demonstrate  a  class  of  feedback  control  laws  which  transform  a  simple 
mechanical  system  with  symmetry  with  exogenous  forces  (controls)  into  a  gyroscopic  system  with 
symmetry,  cf.  Theorem  4.4  below.  Although  the  concept  of  a  Lagrangian  system  with  symmetry 
is  by  now  well-known  (see  [1],  and  [5]),  in  the  interest  of  keeping  our  treatment  self-contained,  we 
give  a  rapid  expose  of  the  basic  ideas  around  the  concept  of  a  gyroscopic  system  with  symmetry. 

DEFINITION  4.1 

A  Gyroscopic  System  with  Symmetry  is  a  5-tuple,  ( Q ,  K,  Y,  V,  G) ,  where 

(1)  (Q,  K)  is  a  riemannian  manifold. 

(2)  Y  is  a  vector  field  on  Q  ,  called  a  gyroscopic  field. 

(3)  V  is  a  function  on  Q  ,  a  potential. 

(4)  G  is  a  Lie  group  with  an  action  $  :  G  x  Q  — ^ ►  Q ,  which  leaves  K ,  Y ,  V  invariant  and  is 
referred  to  as  the  symmetry  group. 

(5)  The  associated  Lagrangian  L  :TQ  It  is  given  by 

L(vx)  =  \K(x){vx,  vx)  +  I((x)(vx,  Y(x))  -  V(x).  (4.1) 


a 

On  the  other  hand,  in  the  framework  of  hamiltonian  mechanics,  a  gyroscopic  system  with 
symmetry  is  characterized  by  a  Hamiltonian  H  :  T*Q  — »•  1R  in  the  following  form, 

H(ax )  =  ^  <  a*  -  K\Y(x)),  ax  -  K\Y(x))  >t-q  +  V(x).  (4.2) 

where  <  •,  •  >t*q  is  the  induced  metric  on  T*Q  defined  in  (2.2). 

The  word  “gyroscopic”  comes  from  the  second  term  in  the  Lagrangian  (4.1),  which  includes 
the  gyroscopic  field  Y .  This  term  is  linear  in  the  velocity  variables  and  is  responsible  for  the 
paradoxical  behavior  of  gyroscopes.  The  Coriolis  force  in  a  rotating  reference  system  and  the 
magnetic  force  due  to  electric  currents  are  examples  of  the  effect  of  gyroscopic  terms  in  the 
Lagrangian  function.  To  see  how  the  gyroscopic  term  enters  the  dynamical  equations,  we  restrict 
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our  attention  for  the  moment  to  a  gyroscopic  system  (without  symmetry  consideration)  on  Hn  (or 
in  local  coordinates)  described  by  the  Lagrangian 


L( x,  v)  =  —  <  v,  M(x)v  >  +  <  Y(x),  v  >  —  V(x), 


(4.3) 


where  M(x )  is  a  symmetric  positive-definite  second-order  tensor,  vectors  x,  v(=  x)  are  in  Hn  ,  Y 
is  a  map  from  lRn  to  ]Rn,  and  V  is  a  real-valued  function.  The  notation  <  -,  •  >  denotes  the  inner 
product  on  the  Euclidean  space  IR.n .  This  is  a  gyroscopic  system  in  the  sense  of  Definition  4.1  with 

K(x)(v,  v )  =  vTM(x)v ,  and  Y(x)  =  M(x)~lY( x). 


Abstractly,  Y  should  be  regarded  as  a  1-form  in  T*Q . 

To  obtain  the  dynamical  equations  associated  with  the  Lagrangian  in  (4.3),  we  invoke  the 
classical  Euler-Lagrange  equations,  cf.  (3.13).  First,  we  find 


ax 

dv 


=  M(x)  ■  v  +  Y(x). 


By  taking  time  derivatives,  we  get 

d  8L 


„r,  N  .  (dM.  ,  \  dY  ,  . 

dt  dv  =  +  (■ a7(l) ' v)  ' v  + 

where  dM/dx  is  a  third-order  tensor,  and  dY/dx  is  a  second-order  tensor.  With  standard 
notations  in  Tensor  Algebra,  cf.  e.g.  [3],  the  above  equation  can  be  rewritten  as 


ddL  .  .  dM.  .  ,  dYf  , 

__  =  #(,).,  + 

Now  we  compute  the  partial  derivative  of  L  with  respect  to  x.  By  definition, 


d_ 

dx 


(vT  M  (x)v) 


w  = 


d_ 

de 


v1'  M  (x  +  ew)u, 


e=0 


d 

de 


M(x  -few):  (uv) 


6=0 

dM 


dx 


( x ) • w  :  (vv)  = 


dM 

dx 


(*)* :  ( vv ) 


w. 


where  (dM  I  dx)*  is  the  cyclic  transpose  of  dM/dx  defined  through 

dM  .  a  dM  ^ 

-(x)  -u-v-w  =  — — (x)-w-u-v,  V  u.v,w  6  M 
v  dx 


dx 


Accordingly,  we  have 


d,T 


dM 


gx(v  U(x)v)  -  9i 


(*)*  :  ( VV )■ 


(4.4) 


(4.5) 


On  the  other  hand, 
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■  w 


d_ 

dx 


( vtY (a:))  •  w 


d_ 

de 


6=0 

dY_ 
dx 


vtY( x  jew) 


(x)  •  W  ■  V  — 


d_Y_ 

dx 


(. x)T  ■  v 


where  the  superscript  T  denotes  the  transpose  of  the  second-order  tensor.  Thus  we  have 


d 


{vT?{x)) 


dr 


dxyu  *  v~"  ~  dx 


( x)T  ■  V. 


(4.6) 


By  substituting  (4.4),  (4.5),  and  (4.6)  in  the  Euler- Lagrange  equations,  we  obtai 


am 


8M 


8Y 


1  dM 


dY 


dV 


M(x)-i  +  ^.):W  +  ^(x)-v  =  (x)‘:(H  +  ~(xf  ■  v  -  ~(x). 


which  can  be  further  written  as 
, ,,  .  .  (1  dM 


dM 


dY  ,  „  dY  ,  ,T 


dx(x)  :{VV)  -  a) 


’  -  ^(X)' 


Define 


(4.7a) 

(4-76) 


The  equations  of  motion  for  a  gyroscopic  system  with  Lagrangian  of  the  form  (4.1)  can  be  then 
expressed  as, 

dV 

M(x)-x  =  T-x-x  —  11 -x  -  -q— (x).  (4.8) 

Note  that  71  is  a  skew-symmetric  tensor,  thus  the  second  term  on  the  right  of  (4.8)  gives  the 
gyroscopic  force  in  the  dynamical  equations  as  discussed  in  [nj.  We  remark  here  that  the  component 
form  of  M(x)~lT  is  nothing  but  the  Christoffel  symbol  associated  with  the  geodesic  flow.  Cf.  e.g. 
[i].  By  multiplying  both  sides  of  (4.8)  with  M~l ,  we  get 

rv-r/ 

x  —  M(x)~1T-x-x  =  M(x)~1lZ-x  —  M(x)~l— — (x).  (4.9) 

( J  X 

In  terms  of  covariant  differentiation,  we  can  write  the  left  hand  side  of  (4.9)  as  V^x.  Moreover, 
with  Y  being  a  1-form,  by  taking  exterior  derivative  of  Y ,  we  obtain,  in  local  coordinates, 
dY(x)(v,  w)  =  ( DY(x)'v)-w  —  ( DY(x)-  w)-v 

=  DY(x)‘V-w  -  DY[x)t-v-w  —  IZ-v-w 

Thus  we  may  write 

1Z  ■  x  =  dY(x)(x ,  ■), 
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which  is  actually  a  1-form.  On  the  other  hand,  the  third  term  on  the  right  of  (4.8)  also  corresponds 
to  a  1-form,  namely,  dV .  Recall  that  K 8  transforms  a  1-form  to  a  vector  field.  We  define  the 
following  notations 

grady  =  A'8  •  dV, 

(curviy)“  =  A8  •  dY( x)(x,  ■). 

Here  grady  denotes  the  gradient  of  V .  As  a  consequence,  (4.9)  can  be  expressed  in  its  invariant 
form  as  follows. 

THEOREM  4.2 

The  invariant  form  of  the  equations  of  motion  of  a  gyroscopic  system  ( Q,K ,  Y,  V )  is 

Vi  a:  +  (curv^Y)8  =  -gradY.  (4.10) 


1 

These  are  the  equations  of  motion  of  a  charged  particle  in  a  magnetic  field  dY .  They  are  thus  a 
special  (simplest,  abelian)  case  of  Wong’s  equation  for  the  motion  of  a  particle  in  a  Yang-Mills  field 
[53],  [17].  The  term  (curViY)8  is  the  corresponding  force  of  interaction  with  the  field. 

EXAMPLE  4.3 

We  consider  the  dynamical  system  treated  in  [11]  in  the  following  form. 

x  =  - ax  -  gy,  y  =  -(3y  +  gx.  (4.11) 


The  skew  terms  in  velocities  —  gy  and  gx  constitute  the  gyroscopic  forces  which  do  no  net  work 
but  affect  the  stability  of  the  system.  It  is  easily  checked  that  this  is  a  gyroscopic  system  with  the 
Lagrangian  in  the  form  (4.3)  with  the  following  entities, 

?)•  (0 1)  (; 


M(x,y)  = 


1  0 
0  1 


Y{x,y)  = 


i 

With  the  equations  of  motion  of  a  gyroscopic  system  (4.8),  or  (4.10),  we  are  ready  to  state 
the  main  ideas  of  gyroscopic  control.  This  notion  is  isolated  here  to  highlight  the  role  of  the 
gyroscopic  term  from  the  viewpoint  of  designing  control  algorithms.  A  simple  mechanical  system 
with  symmetry  with  exterior  forces  can  be  transformed  into  a  gyroscopic  system  with  symmetry 
by  using  suitable  feedback  laws,  that  we  refer  to  as  gyroscopic  feedback.  This  process  is  described 
in  the  following  theorem. 
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THEOREM  4.4 


Consider  a  simple  mechanical  system  with  symmetry,  ( Q,K,V,G ).  Let  the  exterior  force 
exerted  on  this  system  be  denoted  by  a  horizontal  1-form  (a,  0).  Let  Y  be  any  G -invariant 
1-form  on  Q  .  Then,  with  the  feedback  law  ( gyroscopic  feedback ), 

a(vx)  =  -dY(vx,-),  (4.12) 

the  closed-loop  system  becomes  a  gyroscopic  system  with  symmetry  ( Q,K,Y,V,G )  where, 

Y  =  Xs -Y.  (4.13) 

Proof 

We  prove  this  theorem  in  local  coordinates  where  the  riemannian  metric  is  expressed  as 

K(x)(v,  to)  =  vT  M{x)w. 


Also  the  feedback  law  (4.12)  can  be  written  in  local  coordinates  as 

f  d  Y  dYT  \ 

v)  =  ~  (  to (x)  -  {x))  ”•  (4’14) 

Recalling  the  derivation  of  equation  (4.8),  the  dynamical  equations  for  (Q,  K,V,G)  with  exterior 
force  can  be  found  to  be,  in  local  coordinates, 


M(x)  ■  x 


T  ■  x  ■  x 


dV,  s 
to (x)  +  “• 


(4.15) 


where  T  is  defined  in  (4.7b).  With  the  feedback  law  (4.14),  it  is  then  easy  to  see  that  (4.15) 
becomes  (4.8)  which,  in  turn,  corresponds  to  a  system  with  Lagrangian  in  the  form  of  (4.3).  With 
the  transformation  rule  (4.13)  expressed  in  local  coordinates, 

y(x)  =  M(x)-1  Y(x), 


the  system  can  be  further  identified  as  a  gyroscopic  system  with  symmetry,  (Q,  K,Y,  V,  G) .  The 
G -invariance  property  of  Y  follows  from  the  G -invariance  of  the  riemannian  metric  and  the  1- 
form  Y . 

1 

Accordingly,  we  have  a  family  of  gyroscopic  feedback  laws  induced  by  G -invariant  1-forms. 
The  techniques  used  for  analyzing  gyroscopic  systems  with  symmetry  can  then  be  applied  to  study 
the  corresponding  closed-loop  system.  In  particular,  the  method  for  stability  analysis  based  on  the 
energy-momentum  method  which  will  be  developed  in  the  following  is  applicable.  The  gyroscopic 
term  affects  the  dynamical  behavior  in  many  ways.  For  example,  it  changes  the  location  of  equilibria 
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as  well  as  their  stability  properties.  As  a  consequence,  suitable  gyroscopic  feedbacks  may  be  chosen 
to  fulfill  design  objectives.  Much  work  remains  to  be  done  on  general  methods  for  selecting  Y .  The 
dual-spin  problem  illustrated  below  gives  a  simple  example  of  gyroscopic  systems  with  symmetry. 

EXAMPLE  4.5  (Dual-spin  Problem) 

Consider  the  system  consisting  of  a  rigid  body  (platform)  with  on-board  rigid  symmetric 
rotors  moving  in  free  space,  cf.  Figure  4.1.  With  the  rotors  spinning  at  constant  rates  relative  to 
the  platform,  the  dynamical  behavior  can  be  captured  by  a  gyroscopic  system  with  symmetry. 


First  we  consider  the  system  dynamics  with  locked  rotors.  We  assume  that  the  center  of  mass 
of  the  system  is  fixed  in  some  inertial  frame.  Let  B  £  .90(3)  denote  the  orthogonal  transformation 
from  the  body  frame  to  the  spatial  frame,  and  then  describe  the  attitude  of  the  body.  We  have 
B  —  B(l ,  where  the  operator  ~  is  defined  in  (3.14)  and  D  is  the  instantaneous  angular  velocity  of 
the  body  relative  to  the  body  frame.  Let  J  be  the  total  moment  of  inertia  of  the  platform  with 
locked  rotors.  The  Lagrangian  for  this  locked  system  can  be  written  as 

l(b,  bCi)  =  |(n,  jo). 

With  this  Lagrangian,  we  now  invoke  Lagrange-d’Alembert  Principle  specialized  to  the  current 
situation,  cf.  (3.18).  Using  the  parametrizations  introduced  in  Section  3,  the  dynamical  equations 
are  nothing  but  the  Euler’s  equation  for  rigid  body  dynamics, 

j  si  =  -fix  j  a 

Next  we  let  the  rotors  spin  at  constant  rates.  Let  0; ,  i  =  1,2,3  denote  the  relative  angles  between 
the  three  rotors  and  platform,  respectively,  and  0  =  (0i,  92,  #3).  Let  the  corresponding  moments 
of  inertia  of  rotors  relative  to  the  spinning  axis  be  denoted  by  (Is,  )i  1  *  =  1/2, 3  ,  respectively.  The 
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reaction  force  exerted  on  the  platform  from  the  rotors  can  be  derived  from  the  following  gyroscopic 
1-form, 

Y(B)  =  5(1*6),  (4.16) 

where  I5  =  diag{  (Isji,  (Is2)2,  (Is3)3  }•  In  fact,  by  using  the  formula  (4.12),  it  can  be  shown 
that 

a(B,B(l)  =  -5(ftx750). 

With  this  exogenous  force,  the  dynamical  equations  for  the  closed-loop  system  become, 

J!)  =  -fi  X  (Jft  +  I50).  (4.17) 


It  is  readily  checked  that  this  system  is  a  gyroscopic  system  with  symmetry  with  the  following 
entities, 


Q  =  50(3),  K(Buu  Bii2)  =  («l7  Ju2>, 

Y(B)  =  B(  J-”iF0),  V(B)  =  0,  G  =  50(3) 


(4.18) 


where  Bu\ ,  Bu^  £  TbS0( 3).  The  group  action  here  is  G  X  Q  -+  Q ,  (R,  B )  RB,  and  the 

Lagrangian  is 

L(B,  BQ)  =  Jfi)  +  (ft,  Is0).  (4.19) 


| 

Keeping  in  focus  our  program  of  understanding  the  closed-loop  behavior  of  a  system  with 
gyroscopic  feedback,  we  return  to  the  abstract  framework  of  gyroscopic  systems  with  symmetry.  We 
remark  first  that  a  simple  mechanical  system  with  symmetry  in  the  sense  of  Smale  [42]  is  a  special 
case  of  a  gyroscopic  system  with  symmetry.  We  simply  take  Y  —  0  and  consider  the  quadruple 
(Q,  K,  V,  G) .  Many  key  results  in  the  category  of  simple  mechanical  systems  with  symmetry  can  be 
extended  to  gyroscopic  systems  [46].  First,  for  a  gyroscopic  system  with  symmetry,  the  Lagrangian 
(4.1)  is  invariant  under  the  tangent  lift  4>T,  which  follows  directly  from  the  invariance  of  the 
metric  K ,  the  gyroscopic  field  Y ,  and  the  potential  V .  The  Legendre  transformation,  cf.  (3.4),  is 
given  by 

Il(vx)  =  K\vx  +  Y(x)).  (4.20) 

Its  inverse  can  be  then  found  as,  for  ax  6  T*Q , 

(■l1  (ax)  =  Ka(ax)  -  Y(x). 
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It  follows  that  If,  is  a  diffeomorphism,  and  L  is  hyperregular.  As  a  consequence,  the  space 
(TQ,  f II  -  - d&L )  is  a  symplectic  manifold,  where  the  symplectic  form  f lL  is  defined  as  in 
(3.3)  through  the  1-form  ©l,  which  in  turn  can  be  written  as 

•  (w,  ro)  =  K(x)(v  +  Y(x),  u ). 

From  Lemma  4.6,  the  group  G  acts  on  TQ  through  the  tangent  lift  $T  as  a  symmetry  group. 
It  can  be  further  verified  that  this  action  is  symplectic,  namely,  (4>t)*L!l  =  fl:L.  Within  this 
framework,  a  momentum  mapping  J  :  TQ  ->  Q *  can  be  constructed  such  that  the  infinitesimal 
generator  of  the  action  <j?T  corresponding  to  £  6  Q  is  the  vector  field  induced  by  the  function 
(J,  £)  :  TQ  —r  R,  through  the  symplectic  structure,  cf.  (3.9).  Consequently,  we  ha,ve  the 

following  theorem. 

THEOREM  4.6 

The  gyroscopic  system  with  symmetry  ( Q,K,Y,V,G )  has  the  following  properties: 

(i)  The  1-form  corresponding  to  L  defined  in  (3.2)  is  invariant  under  the  tangent  lift,  i.e. 
($J)*  ©L  =  Ql. 

(ii)  There  is  an  associated  Ad* -equi variant  momentum  mapping  J  :  TQ  — ►  Q* , 

J(vx)( 0  =  (lL(vx),  £q(x))x  =  <  vx  +  Y(x),  £q(x)  >x,  (4.21) 

where  £  G  Q  is  an  element  in  the  Lie  algebra  of  G ,  £q(x)  denotes  the  infinitesimal  generator 
of  £  on  Q .  Here  Q*  denotes  the  dual  of  the  Lie  algebra  Q . 

(iii)  The  momentum  mapping  defined  in  (4.21)  is  a  vector-valued  integral  of  any  vector  field 
induced  by  a  G-invariant  function  on  TQ  through  an  analogous  formula  in  (3.9).  In  particular, 
it  is  an  integral  of  the  Lagrangian  vector  field  X hl  ■ 

Proof 

For  (i),  we  note  that  ($J)*T  =  L ,  by  Lemma  4.6.  Since  the  exterior  differentiation  commutes 
with  the  pull-back  operator,  (i)  follows  immediately.  Statements  (ii),  and  (iii)  can  be  shown  by 
directly  applying  Theorem  4.2.2  and  Corollary  4.2.14.  in  [i]. 

I 

The  quadruple  (TQ,  $t,  J )  is  an  example  of  a  hamiltonian  G-space.  The  energy  function 

for  the  gyroscopic  system  can  be  derived  as.  cf.  (3.7), 

r)  —  (^l(w)j  ^ x)x  h  ( Vx ) , 

=  ^  vx  T  I  (x'),vx  x  —  —  rq,  vx  <C  tq.I  (a;)  4-  I  (x), 

=  ~  <  vx,  vx  »*  +  V(x).  (4.22) 
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It  is  easy  to  see  that  the  energy  function  is  not  affected  by  the  presence  of  the  gyroscopic  field  Y . 
However,  the  dynamics  are  different  from  what  one  would  see  if  Y  =  0 .  The  differences  in 
the  dynamical  behavior  arise  from  the  Y -dependent  symplectic  "2-form  fin  -  In  particular,  the 
gyroscopic  term  in  the  Lagrangian  gives  rise  to  the  magnetic  terms  in  the  symplectic  2-form.  On 
the  other  hand,  on  the  momentum  phase  space  T*Q ,  the  Hamiltonian  associated  to  the  system  is, 
cf.  (4.2), 


H(ax)  =  El  o  CL\ax\ 

=  \  <  KHax)  -  Y(x),  K\ax )  -  Y{x)  +  V(x), 

=  \  <  <*x  -  K\Y{x)),  ax  -  IC\Y{x ))  >t-q  +  V(x). 


Accordingly,  on  the  momentum  phase  space,  the  Hamiltonian  is  affected  by  the  gyroscopic  term 
through  a  momentum  shift,  while  the  canonical  2-form  wo  is  unchanged.  This  subtlety  is  best 
explained  by  the  following  example. 


EXAMPLE  4.7 

We  consider  again  the  system  in  Example  4.3.  The  energy  associated  with  (4.11)  on  TQ  is 
HL(x,y,x,y)  =  ~  (x2  +  y2  +  ax 2  +  (iy2)  , 


with  the  symplectic  2-form  in  matrix  representation 

0  g  1  0\ 

[  -9  0  01 


[fin] 


-10  0  0 
0  -1  0  o/ 


This  can  be  checked  from  the  differential  equation  (4.11),  cf.  (3.9b), 


k  =  Ahl(x)  =  ([fiL]-1)2,  VHl, 


where  x  =  (x,  y ,  x,  y) .  The  left-upper  2x2  block  in  [fin]  is  called  the  magnetic  part.  On  the  other 
hand,  on  T*Q ,  we  have  the  conjugate  momentum  variables  defined  by 

Pi  =  x  +  gy,  p2  =  y. 

The  dynamical  equation  (4.11)  can  be  written  as 

x  -  pi  -  gy,  y  =  p-2, 

Pi  =  -  ax,  go  =  -fiy  +  g(pi  -  gy), 

which  is  a  hamiltonian  system  with  the  Hamiltonian  function 
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H(x,y,Pi,P2)  =  x  {{ih-gy)2  +  pi  +  CXX2  +  (3y2)  . 


The  symplectic  structure  is  the  canonical  symplectic  2-form  u>0,  i.e.  in  matrix  representation, 


M  = 


/  o  o  1  0\ 
[  0  0  0  1 
-1  0  0  0 
Vo  -loo/ 


1 

In  summary,  the  gyroscopic  term  affects  the  symplectic  2-form  on  TQ  side  while,  on  T*Q 
side,  it  affects  the  Hamiltonian  function.  To  gain  more  insight  on  how  the  gyroscopic  field  enters 
the  symplectic  structure  Q,l,  we  consider  an  even  simpler  case  than  (4.3).  We  assume  that  the 
second-order  tensor  M  is  independent  of  x  in  (4.3).  It  can  be  easily  seen  that  the  symplectic 
2-form  is  now,  cf.  (3.3), 


fi£(?>v)((«i>  wi),  (u2,  w2)) 

OY  ,  BY 

=  — -  •  Uo  ■  U\  +  M  •  Wo  ■  m  —  — —  •  U\  •  Uo 

ox  ox 


=  («1  Wl) 


M  ■  W\  ■  uo 


The  block  ^  is  the  so-called  magnetic  term. 


5.  Reduction,  Relative  Equilibria  and  Stability 

By  recognizing  the  symmetry,  under  suitable  regularity  hypotheses,  it  is  possible  to  reduce  a 
gyroscopic  system  with  symmetry  ( Q,K ,  Y,  V ,  G)  to  a  lower-order  dynamical  system.  The  reduction 
process  has  a  long  history.  For  Jacobi  and  Liouville[i]  [5],  this  meant  reduction  of  the  Hamilton’s 
equation  via  first  integrals  in  involution.  For  Routh  [37],  this  meant  a  process  of  eliminating 
ignorable  variables.  In  the  following,  we  shall  discuss  the  reduction  from  two  modern  points  of 
view,  namely,  symplectic  reduction  and  Poisson  reduction. 

First,  we  consider  symplectic  reduction  in  the  sense  of  [32].  As  discussed  in  Section  3,  if  L  is 
regular,  (TQ ,  fix,)  is  a  well-defined  symplectic  manifold.  By  the  Property  (i)  in  Theorem  4.6,  the 
Lie  group  G  acts  symplectically  on  (TQ  ,  ).  Also,  from  Property  (ii )  in  Theorem  4.6,  there  is  an 

Acl* -equivariant  momentum  mapping  .J  for  this  action.  Thus  all  the  conditions  in  the  Symplectic 
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Reduction  Theorem,  see  Theorems  4.3.1,  4.3.5,  pp.  299,  304  in  [l]  are  satisfied,  and  we  can  state 
the  following  reduction  theorem  specialized  to  gyroscopic  systems  with  symmetry. 

THEOREM  5.1  (Marsden- Weinstein) 

Consider  the  gyroscopic  system  with  symmetry  (Q,K,Y,V,G).  Assume  that  p  G  Q*  is  a 
regular  value  of  the  momentum  mapping  .J ,  as  defined  in  (4.21),  and  that  the  isotropy  subgroup 
Gy ,  defined  by  Gy  =  {  g  G  G  :  Ad*_1p,  =  p},  under  the  Ad*  action  on  Q*  acts  freely 
and  properly  on  J~1(p),  then  ( TQ)y  =  J~1{p)/Gy,  has  a  unique  symplectic  form  £ly  with 
the  property  7r*  f ly  =  i*  where  rr y  :  J~1(p)  —>  ( TQ)y  is  the  canonical  projection  and 

iy  :  J~1(p)  '-»■  TQ  is  the  inclusion  map.  Letting  Hl  be  as  in  (4.22),  the  flow  Ft  of  Xhl  induces 
a  flow  Ft  on  ( TQ)y  satisfying  Ky-Ft  —  Ff  ■  i rM.  This  flow  is  a  Hamiltonian  flow  on  ( TQ)y  with 
a  Hamiltonian  function  satisfying  H£  -liy  =  Hl  • iy ,  with  respect  to  the  symplectic  structure 

JV 

1 

The  function  H£  on  the  reduced  space  is  called  the  reduced  energy.  The  corresponding 
vector  field  Xjj»  on  the  reduced  space  (TQ)y  is  called  the  reduced  dynamics.  Thus  in  symplectic 
reduction,  we  first  restrict  the  dynamics  to  a  level  set  of  the  momentum  mapping,  and  then  factor 
out  the  isotropy  subgroup. 

Next,  we  consider  Poisson  reduction  [31].  We  first  recall  the  basic  setup  of  Poisson  manifolds. 
A  Poisson  manifold  P  is  a  smooth  manifold  equipped  with  an  IR-bilinear  map  (Poisson  structure) 
on  the  space  of  smooth  functions,  {-,-}p  :  C°°(P)  x  Cco(P) C<X,(P)  satisfying  the  axioms, 
for  f,g£C°°(P), 

(i)  {f,g}p  =  -  {§J}p 

(ii)  {fg-h}p  =  g{f,h}p  +  f{g,h}P 

(iii)  {/>  {g-,h}p}P  +  { g,{h,f}p}p  +  {h,{f,g}p}p  =0. 

Associated  to  a  Poisson  structure,  there  is  a  unique  twice  contravariaj.it  skew- symmetric, 
smooth  tensor  field  A  on  P  such  that  {f,g}P  =  A (df,dg),  where  df ,  dg  are  differentials  of  /, 
g ,  respectively.  The  tensor  field  A  defines  a  vector-bundle  morphism,  A#  :  T*P  — ►  TP ;  ax  ^ 
A*(ax)  G  TXP ,  satisfying, 

Px  (A*(ax))  =  A  (x)  {ax,px)  for  all  (5X  G  T*XP. 

Let  G  be  a  Lie  group  and  let  $  :  Gx  P  — *-  P ,  (g,x)  i->  4rs(a-),  be  a  group  action  such  that 
'$'(,(•)  is  a  Poisson  morphism  for  every  g  G  G,  i.e.  If g  :  P  — *■  P  is  an  isomorphism  and  preserves 
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the  Poisson  structure.  Suppose  that  the  action  is  proper  and  free.  Then  the  quotient  space  P/G 
is  a  manifold  which  carries  a  Poisson  structure  {•,  -}p/g  induced  from  the  one  on  P  satisfying, 
for  f,g€  C°°(P/G) , 

{f,g}p/G°K  =  {/0  7T,  g°ir}p.  (5.1) 

Here  it  :  P  — *  P/G  is  the  canonical  projection.  By  construction,  it  is  a  Poisson  morphism. 

G-equivariant  dynamics  on  P  induce  dynamics  on  P/G.  Suppose  h  :  P  — >  HR,  is  a  G- 
invariant  Hamiltonian  function  on  P,i.e.,  h  (\I/3(x))  =  h(x),  V  g  £  G.  Define  a  vector  field  X h 
through 

xh[f]  =  {/,  h}p ,  V  /  e  G°°(P).  (5.2) 

where  Xh[f]  denotes  the  Lie-derivative  of  the  vector  field  Xh  on  the  function  /.  The  Hamiltonian 
h  descends  to  h:P/G-*l R  and  determines  a  Poisson-reduced  dynamics  X-h  on  P/G  by 

Xh  [/]  =  {/,  h}p,a ,  V/  €  G°°  (P/G).  (5.3) 

Here  h  ([a;])  =  h(x)  for  an  equivalence  class  [x]  in  P. 

Recall  that  the  symplectic  manifold  (TQ,Ql)  has  a  Poisson  structure  induced  from  the 
symplectic  structure,  namely,  for  /,  g  €  C°°(TQ), 

{/,  g}L(vx)  =  df(vx)-Xa(vx)  =  SlL(vx)(TLL(df),  n L(dg)),  (5.4) 

cf.  (3.5).  Since  the  energy  function  Hl  is  G-invariant,  we  carry  out  the  Poisson  reduction  as 
follows.  Assume  G  acts  on  TQ  freely  and  properly.  Let  t  be  the  projection  from.  TQ  to  TQ/G, 
f,g  £  C°°(TQ/G),  the  induced  Poisson  bracket  of  /  and  g  is  defined  analogous  to  (5.1)  as 

{/,  =  {/of,  jof}L.  (5.5) 

Referring  to  the  framework  of  Poisson  reduction,  we  can  identify  the  induced  Hamiltonian  FIl  and 
associated  dynamics  XfjL  as: 

Hl  o  f(vx)  =  Hl{vx ),  (5-6) 

XkL[f]  -  {/,  Hl}>,  V  /  G  C°°(TQ/G).  (5.7) 

Here  the  vector  field  X^  is  called  the  projected  Hamiltonian  vector  field  on  TQ/G. 

The  reductions  discussed  here  are  on  the  lagrangian  side,  or  TQ  side.  We  could  perform  a 
similar  reduction  process  on  T*Q  side,  or  hamiltonian  side,  by  noting  that  the  Hamiltonian  function 
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on  T*Q ,  namely  II  in  (4.2),  is  invariant  under  the  cotangent  lift  $T*  (this  follows  from  (2.6)). 
The  underlying  symplectic  manifold  is  (T*Q,  <n0),  with  the  corresponding  momentum  mapping, 

J  :  T*Q  -»•  6?*;  (J(ax),  $  =  (ax,  £<,(*))*,  V(€5.  (5.8) 

Since  L  is  hyperregular,  the  reductions  on  TQ  and  T*Q  are  equivalent,  but  we  shall  use  the  one 
on  TQ  side  in  the  following  development,  bearing  in  mind  that  the  Lagrange-d'Alembert  Principle 
is  formulated  there. 

We  proceed  to  discuss  a  characterization  of  relative  equilibria.  The  concept  of  relative 
equilibrium  goes  back  to  Poincare.  In  the  context  of  symplectic  reduction,  we  define  the  notion  of 
relative  equilibrium  as  follows,  cf.  Theorem  5.1. 

DEFINITION  5.2 

A  point  vx  in  TQ  is  called  a  relative  equilibrium  if  E  (TQ)m  is  a  fixed  point  for  the 

symplectic-reduced  vector  field  Xjj*  ,  where  p,  =  J(vx). 

I 

In  the  context  of  Poisson  reduction,  we  can  define  a  similar  notion,  cf.  (5.6),  (5.7).  A 
point  vx  in  TQ  is  called  a  relative  equilibrium  for  Xjjl  with  respect  to  the  Poisson  reduction 
if  Xfj  (f(vx))  —  0.  It  turns  out  that  the  two  notions  of  relative  equilibrium  are  equivalent.  In 
fact,  it  can  be  shown,  cf.  [i],  that,  for  both  cases,  vx  is  a  relative  equilibrium  iff  there  exists  a 
£  E  Q  such  that  the  flow  of  Xhl  , 

Fx :H[  (Vx)  =  ^exp(t^)  (w)>  (5-9) 

namely,  the  dynamical  orbit  is  simply  a  group  orbit.  Thus  if  the  observer  were  to  be  set  in  uniform 
motion  according  to  the  one-parameter  group  exp(t£) ,  then  for  such  a  moving  observer,  a  relative 
equilibrium  will  appear  to  be  stationary.  For  instance,  if  G  =  50(3),  then  .Fy  (nr)  corresponds 
to  a  uniform  rotation  about  a  fixed  axis  £  in  space  with  the  rotational  speed  |£| .  In  a  central  force 
field,  a  relative  equilibrium  for  the  motion  of  a  point  mass  corresponds  to  a  circular  orbit,  cf.  [49]. 
Relative  equilibria  can  be  characterized  by  the  following  result  of  Souriau-Smale-Robbin. 

THEOREM  5.3 

vx  G  TQ  is  a  relative  equilibrium  for  X 'h  iff  there  exists  a  f  E  Q  such  that  vx  is  a  critical  point 
of  £),  where  (J,  £):TQ  — >  R  is  the  real- valued  function  given  by  vx  (J{vx),  £), 

associated  to  the  the  momentum  mapping  J . 

H 

In  particular,  for  gyroscopic  systems  with  symmetry,  we  have,  cf.  (5.8),  (4.22), 
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(5.10) 


Ht(vx)  =  |  <  vx,  vx  >x  +  V(z)-  <  vx  +  Y(x),  £Q(x)  >x, 

=  ^  <  t?*  -  £q(x),  -  fg(a:)  »* 

+  1  (a;)  -  <C  1  (a;),  £<3(2;)  "^>x  ~  ^  “C  £q(x)i  £q{x) 

From  Theorem  5.3,  it  is  then  easy  to  check  that  the  necessary  and  sufficient  conditions  for  vx  to 
be  a  relative  equilibrium  are 

vx  -  (q(x),  (5.11) 

and 

dx[V(x)  —  <  y(a;),  £q(x)  —  —  <C  £q{x)-  £q(x)  —  0. 

We  thus  have  the  following  algorithm  ( principle  of  symmetric  criticality  [35])  to  find  relative 
equilibria. 

ALGORITHM  5.4 
0.  Pick  teg. 

1.  Search  for  the  critical  points  xe  of  the  function 

V^-.Q  ->■  R;  V^x)  =  V(x)  -  <  Y(»),  £q(x)  >x  -  l-  <  $q(x),£q(x)  >x  .  (5.12) 

2.  Substitute  xe  in  (5.11)  to  find  the  corresponding  ve  =  £g(a:e). 

I 

We  note  that  the  computation  in  step  1  is  fully  on  the  configuration  space.  Thus  the  process 
of  searching  for  a  relative  equilibrium  is  greatly  simplified.  We  remark  that,  for  simple  mechanical 
systems  with  symmetry,  the  principle  of  symmetric  criticality  stated  above  appears  as  Theorem  1.1 
in  Part  II  of  Smale[42].  Smale  also  notes  that  special  versions  have  been  known  earlier,  e.g.  in 
the  study  of  symmetric  geodesics.  See  also  p.  355  of  [1],  Theorem  16.7  in  Hermannfig],  Arnold[6], 
and  Palais[35],  Here  the  augmented  potential  function  V?  has  one  additional  term  to  accommodate 
the  gyroscopic  effects.  Through  this  term,  we  can  change  the  number  and  locations  of  the  critical 
points.  This  provides  us  an  effective  tool  in  controlling  the  phase  portrait.  Cf.  [46]. 

There  is  an  additional  symmetry  in  the  augmented  potential  \\  .  First,  we  define  the  stabilizer 
of  f  to  be  G $  =  {g  £  G  |  Adg(l ;)  =  £}  C  G,  where  Ad  is  the  adjoint  action  of  G  on  Q  defined 
in  (2.3).  G £  is  actually  a  subgroup  of  G ,  and  thus  defines  an  action  on  Q  .  By  an  argument  similar 
to  the  one  in  the  proof  of  Lemma  4.6,  it  can  be  shown  that  is  invariant  under  the  action  of  G ^ 
on  Q ,  i.e. 


26 


V^g(x))  =  Vt(x),  VgeGt. 


(5.13) 


We  assume  that  the  quotient  space  Q/G^  is  well  defined.  Denote  the  projection  from  Q  to  Q/G % 
by  7 .  By  (5.13),  we  can  define  an  induced  function  on  Q/G £  from  the  augmented  potential 
such  that  the  diagram  in  Figure  5.1  commutes,  namely, 

V*  o  ■Ki  =  Vv  (5.14) 


Figure  5.1.  Symmetry  of 


This  symmetry  of  Vj  will  be  used  later  in  establishing  a  stability  result  associated  with  . 

As  mentioned  earlier,  reductions  could  be  worked  out  on  the  Hamiltonian  side  as  well.  Thus 
there  is  a  similar  algorithm  corresponding  to  Algorithm  5.4  on  the  T*Q  side.  We  only  need  to 
find  the  corresponding  conjugate  momentum  variable  pe ,  by  substituting  xe  obtained  in  Step  1  of 
Algorithm  5.4  in  the  formula,  pe  =  Kb  (Y( xe)  —  ^g(se))  •  The  point  (xe,  pe)  in  the  momentum 
phase  space  T*Q  is  then  a  relative  equilibrium  corresponding  to  the  reduction  on  T*Q  with  respect 
to  the  cotangent  lift  action. 

We  now  address  the  stability  of  relative  equilibria.  Although  both  symplectic  reduction  and 
Poisson  reduction  lead  to  equivalent  notions  of  relative  equilibria,  the  associated  stability  issues  are 
quite  different.  In  the  following,  we  shall  state  the  main  ideas  in  a  somewhat  more  general  setting 
than  needed  for  this  paper.  However,  we  think  that  this  added  generality  keeps  the  treatment  free  of 
confusing  details.  In  general,  let  B,  P  be  differentiable  manifolds,  and  G  be  a  Lie  group.  Consider 
a,  principal  G -bundle,  (P,  G,  B) ,  namely,  G  acts  differentiably  on  P  freely  and  properly,  B  ==  P/G 
is  the  quotient  space  of  P  with  the  canonical  projection  7r  :  P  — ►  B  being  differentiable.  Moreover, 
P  is  locally  trivial,  that  is,  every  point  u  £  B  has  a  neighborhood  U  such  that  there  is  a  mapping 
from  7r-1(C/)  to  U  X  G ,  z  (t t(z),  <b{z))  which  is  a  diffeomorphism  and  o[g  ■  z)  =  g  ■  <f>(z) ,  for 
all  g  6  G .  See  Figure  5.2  for  an  illustration  of  the  geometric  structure  of  such  an  object.  For  more 
details,  see,  e.g.  [34]. 
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Figure  5.2.  Principal  G -bundle 


A  vector  field  X  on  P  is  said  to  be  projectable  if  for  each  /  6  P(B) ,  there  exists  a  /  £  P(B) 
such  that  =  cf,  e.g.  (30],  [18].  Now.  given  a  projectable  vector  field  X  on  P, 

the  corresponding  projected  vector  field  X  on  B  is  defined  in  the  following  way.  Given  a  smooth 
function  /  on  B ,  the  Lie  derivative  of  X  on  f  is  defined  as 

X[  /  1  =  /,  or  X]  /  ]  =  v  =  X[  fox].  (5'15) 

It  is  easy  to  verify  that  the  vector  field  Xe  defined  in  (5.2)  is  projectable  with  the  projected  vector 
field  X-h  defined  in  (5.3)  in  the  above  sense. 

DEFINITION  5.5 

For  the  principal  G-bundle,  (P.G.B),  a  point  zEfis  called  a  relative  eyuilibnnm  of  a 
projectable  vector  field  lEl(P)  if  zr(s)  is  an  equilibrium  of  the  associated  projected  vector  field 
Moreover,  a  relative  equilibrium  zEfis  relatively  stable  module  G  if  the  equilibrium 
Ti{z)  is  Lyapunov  stable  with  respect  to  the  projected  vector  field  X . 


REMARK  5.6  ,  ,  . 

m  [28],  the  smooth  manifold  structure  of  the  quotient  space  P/G  is  not  explicitly  invoked  m 

defining  the  notion  of  stationary  motion  and  relative  stability  modulo  G.  However,  when  the  group 

action  is  free  and  proper,  P/G  is  a  manifold.  This  is  the  case  considered  in  this  paper,  and  hence 


Definition  8.13,  p.  242  in  [28]  is  equivalent  to  Definition  5.5. 

■ 

For  a  gyroscopic  system  with  symmetry,  the  definition  of  relative  equilibrium  vx  €  TQ 
in  Definition  5.4  matches  with  the  Definition  5.5  by  noting  that  the  principal  G-bundle  is 
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now  (TQ,G,TQ /G) .  Accordingly,  the  relative  equilibrium  vx  is  relatively  stable  modulo  G 
in  TQ  if  f{yx )  is  a  stable  equilibrium  with  respect  to  the  projected  Hamiltonian  vector  field 
Xfj  .  On  the  other  hand,  in  the  symplectic  reduction  process,  we  have  the  bundle  structure 
(J~1(p),  G^,  (TQ)fj.).  The  relative  equilibrium  defined  in  Definition  5.2  can  be  regarded  as  a 
relative  equilibrium  with  respect  to  this  principal  G -bundle.  Correspondingly,  we  may  define 
relative  stability  modulo  G ^  in  with  respect  to  the  reduced  dynamics  XH*  .  Since  the 

space  (TQ) fj,  is  diffeomorphic  to  a  symplectic  leaf  in  TQ/G ,  relative  stability  modulo  G  in  TQ 
implies  relative  stability  modulo  G^  in  J~1(p).  The  converse  is  illustrated  by  the  following  theorem 
from  [28],  Theorem  8.17,  p.  244,  see  also  [50]  [24]. 

THEOREM  5.7 

Let  v ®  be  a  relative  equilibrium,  cf.  Definition  5.2.  Definiteness  of  the  Hessian  D2  H £  at 
Tn(vx)  £  (TQ)n  implies  the  relative  stability  modulo  G  in  TQ  of  vx ,  if  there  exists  a  neighborhood 
W  of  f(v%)  G  TQ/G  such  that  the  rank  of  the  Poisson  structure  {•,  •}< ,  defined  in  (5.5),  is  constant 
in  W. 

I 

Those  points  vx  in  TQ  satisfying  the  constant-rank  condition  stated  in  the  above  theorem  will 
be  referred  to  as  generic  points.  The  following  example  demonstrates  that  the  sufficient  condition 
in  Theorem  5.7  is  essential.  This  example  is  from  [28].  A  detailed  discussion  can  be  also  found  in 

[24]. 

EXAMPLE  5.8 

Consider  a  symplectic  manifold  ( P ,  u) ,  where 

P  =  IR4  =  {  (?1,  q2,  Pl,  p2)  },  u  =  dqi  A  dpi  +  dq2  A  dp2.  (5.16) 


Let 

G  =  A//+(IR)  =  {  (a,  6)  G  IR2  }  with  the  group  structure 

(5.17) 

(a,  b )  •  (c,  d )  =  (a  -f  c,  b  +  ea  d). 

It  can  be  shown  that  G  defined  in  (5.17)  is  a  Lie  group.  We  define  the  action  of  G  on  P  as 
G  X  P  -)■  P;  ((a,Z>),  (qi,q2,Pi,P2))  >—  («  +  7i>  b  +  ea  q2,  pi,  e~a  p2). 

It  is  easy  to  check  that  this  is  a  symplectic  action  on  P .  This  action  is  also  free  and  proper.  It 
follows  that  P/G  is  a  manifold  (~  IR2).  The  symplectic  structure  u  in  (5.16)  defines  a  Poisson 
bracket  on  P(P)  which,  in  turn,  induces  a  Poisson  structure  on  P/G .  Let  a  Hamiltonian  function 
H  be  defined  as  H(qi,q2,pi,P2)  =  P2  e91,  which  is  a  G -invariant  function.  It  can  be  checked 
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that  (0,2,  0,0)  E  P  is  a  relative  equilibrium  corresponding  to  the  vector  field  Xh  ■  Moreover,  this 
relative  equilibrium  is  relatively  stable  modulo  G ^  in  J-1(/i),  since  the  quotient  space  J'1(fj,)/Gll 
degenerates  to  a  point.  However,  it  has  been  shown  in  [24]  that  this  relative  equilibrium  is  not 
relatively  stable  modulo  G  in  P .  Note  also  that  that  the  induced  Poisson  structure  doesn’t  have 
a  constant  rank  at  (0,2, 0,0)  and  hence  the  condition  in  Theorem  5.7  does  not  hold. 

I 

There  are  several  methods  for  determining  relative  stability  in  the  appropriate  sense.  For 
example,  the  Energy- Casimir  method  [20]  [46],  or  the  Lagrange-Multiplier  method  [29]  [49]  [46] 
can  be  used  to  determine  relative  stability  modulo  G  in  TQ .  On  the  other  hand,  the  Energy- 
Momentum  method  [39]  [41]  is  useful  in  determining  relative  stability  modulo  G M  in  J_1(/i).  For 
simplicity,  we  will  drop  the  underlying  spaces  in  the  definition  of  relative  stability,  e.g.  we  say 
merely  relative  stability  modulo  G The  underlying  space  is  clear  from  the  context.  In  the  next 
section,  the  Energy-Momentum  method  will  be  adopted  to  study  stability  properties  of  relative 
equilibria  for  gyroscopic  systems  with  symmetry. 


6.  Energy-Momentum  Method  for  Gyroscopic  Systems 

As  pointed  out  in  Section  4,  the  use  of  gyroscopic  feedback  laws  can  affect  the  location  of 
relative  equilibria  and  their  stability  properties.  In  the  work  of  Bloch,  Krishnaprasad,  Marsden,  and 
Alvarez  [7],  an  example  of  rigid  body  stabilization  using  such  a  control  law  is  considered.  Here  we 
give  a  general  method  to  explore  stability  under  gyroscopic  feedback  laws.  A  key  requirement  is  to 
obtain  stability  criteria  that  are  explicit  in  the  parameters  of  the  feedback  lawq  e.g.  the  gyroscopic 
field. 

In  this  section,  the  relative  stability  modulo  will  be  examined  via  the  energy-momentum 
mapping.  Here  we  apply  the  energy- momentum  method  to  the  general  framework  of  gyroscopic 
systems  with  symmetry.  The  block- diagonalization  techniques  for  simple  mechanical  systems  with 
symmetry  is  extended  here  to  account  for  gyroscopic  terms.  The  decomposition  of  the  symplectic 
structure  is  also  presented.  Key  references  for  this  Section  are  [39]  [41]. 

Let  (P,  u))  be  a  symplectic  manifold  on  which  the  Lie  group  G  acts  symplectically  and 
let  J  :  P  —  Q*  be  an  Ad*  -equivariant  momentum  mapping  for  this  action  (see  Section  2  for 
definitions).  Assume  we  could  perform  symplectic  reduction  on  P  in  the  sense  of  Marsden  and 
Weinstein  [32].  The  reduced  phase  space  is  denoted  by  P^  =  J~1(n)/Gfj.  Let  H  :  P  —  1R 
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be  invariant  under  the  action  of  G.  It  induces  a  Hamiltonian  function  H M  on  P M  satisfying 
o  7tm  =  Ho  i where  7rM  :  J_1(/z)  — >  P^  is  the  canonical  projection  and  ^  P  is 

the  inclusion  map.  We  are  interested  in  the  stability  property  of  a  relative  equilibrium  associated 
with  the  reduced  dynamics  on  the  reduced  space  PM,  or  the  relative  stability  modulo 

in  By  construction,  is  a  first  integral  of  the  reduced  dynamics.  Thus  if  has  a 

strict  local  minimum  at  Xn(ze)  where  2e  is  a  relative  equilibrium,  then  jEP  serves  as  a  Lyapunov 
function.  Standard  Lyapunov  stability  analysis  can  be  applied  to  conclude  stability.  Since,  for 

2  €  Jr_10)  C  P, 

ff'W*))  = 

the  condition  for  7rM(2e)  being  a  strict  local  minimum  of  is  equivalent  to  the  condition  for 
ze  being  a  strict  local  minimum  of  H modulo  the  tangent  directions  of  the  group  orbit, 
■  z .  This  in  turn  corresponds  to  checking  that  the  relative  equilibrium  2e  solves  the  constrained 
minimization  problem, 

minimize  H(z)  subject  to  J(z)  —  j.Le  =  J(ze). 

This  problem  could  be  further  formulated  as  checking  ze  to  be  a  strict  local  minimum  of  H  —  { J,  £) 
in  all  directions  on  )  except  along  the  tangent  directions  to  the  group  orbit  generated  by 

G n ,  where  £  serves  as  the  Lagrange  multiplier.  These  heuristic  remarks  are  formalized  in  the 
following,  giving  rise  to  the  energy-momentum  method,  cf.  [40],  [41],  [36],  [39]. 

Define  the  energy-momentum  functional 

H&)  =  H(z)-(J{z),  £).  (6.1) 

From  the  relative  equilibrium  theorem,  cf.  Theorem  5.3,  each  relative  equilibrium  of  the  system  is 
a  critical  point  of  H $  ,  for  some  £  E  Q ,  namely, 

DH^(ze)  ■  6z  =  0,  Vfe  G  TZP. 

From  previous  discussions,  the  definiteness  of  the  second  variation  of  H $  on  a  subspace  S  of  2b.  P 
satisfying 

S  =  TZcJ~l(/.ie)  /  T*4(GM  •  ze),  (6.2) 

implies  the  relative  stability  modulo  GM  of  the  relative  equilibrium  ze .  One  way  to  find  such  a  space 
S  is  to  construct  a  complement  of  TZt{Gil-ze)  in  Tb.  J_1(/re)  such  that  TZc  )  =  S  ®  TZc  ( G ^ 

ze).  Sin ce  =  Iver  D.J(ze),  which  is  the  kernel  of  the  operator  D.J{ze),  we  summarize 

the  energy  momentum  method  for  relative  stability  as  follows. 
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ALGORITHM  6.1  (Energy-Momentum  Method) 

0.  Pick  f  6  Q . 

1.  Solve  the  problem  DH^(z)  ■  Sz  =  0,  V  Sz  E  TZP ,  for  a  relative  equilibrium  ze. 

2.  Compute  fie  —  J(ze)  and  determine  the  space  Ker  DJ(ze). 

3.  Find  S  C  Ker  DJ(ze )  such  that  Ker  DJ(ze )  —  S  ®  TZr  (G ^  ■  ze). 

4.  Check  the  second  variation  of  H ^  on  <S.  Definiteness  of  the  second  variation  implies  stability. 

I 

For  visualizing  the  geometric  picture,  see  Figure  6.1. 


Now  we  restrict  our  consideration  to  gyroscopic  systems  with  symmetry  introduced  in 
Section  4.  The  underlying  space  is  P  —  TQ  with  the  symplectic  structure  fix.  ■  In  this  setting,  the 
momentum  mapping  is  given  by,  cf.  (4.21), 

J{vx)(0  =  <vx+Y(x),  £q(x)>x,  (6.3) 

and  the  energy-momentum  functional  is  JI^(vx)  =  Ii\(vx)  +  V$(x),  where,  cf.  (5.12), 

Kz{vx)  =  ^  <  VX  -  £q(x),  Vx  -  £,q(x)  >x, 

V^x)  =  V(x)~  <  Y(x),  £q(x)  »*  <  £q(x),?q(x)  »r  . 

We  define  the  space 
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(6.4) 


Afx  =  {  i?q(x)  :  r)£G  }, 

which  is  a  subspace  of  TXQ  ,  and  thus  Af  =  UX&Q  Afx,  is  a  subbundle  of  TQ  .  We  then  decompose 
TXQ  =  j\fx  ®Afx  where  Afx  is  the  orthogonal  complement  of  A fx  with  respect  to  the  inner  product 
associated  with  the  riemannian  metric.  Every  element  v  E  TXQ  can  be  thus  written  uniquely  as 
V  =  t)q(x)  +  v,  for  ? 1  £  G,  V  e  Afx  •  With  this  decomposition,  the  function  IQ  could  be  further 
written  as, 

v )  =  \vq(x)  ~  foC*)!!2  +  g  II2- 

Note  that  at  relative  equilibrium  (xe,  £q(xe)),  we  have  ve  —  ^g(xe),  and  1  =  0.  Thus  the  second 
term  in  IQ  is  nonnegative  and  vanishes  at  a  relative  equilibrium  with  a  positive  semi-definite 
second  variation.  Define 

Hi(x,  t?)  =  ~  j|  vq(x)  -  €q(x)  II2  +  vdx)-  (6-5) 

From  the  above  observations,  relative  stability  modulo  G ^  can  be  determined  from  this  modified 
function  on  the  space  of  Q  x  G  •  Moreover,  taking  variations  in  the  space  Q  x  G  corresponds  to 
taking  variations  in  the  subbundle  Af  CTQ.  We  define  the  following  embedding 

2  :  Q  x  Q  —  TQ:  (x,  rj)  >—  (x,  tjq(x)). 


Hence  IQ  =  Qo2.  A  G -action  on  the  space  Q  x  G  can  be  constructed  as 

$  :  Gx{QxG)  —  QxG 

( g ,  (x,  v))  (9-xi  Adg  rj), 

where  Ad  is  the  adjoint  action  defined  in  (2.3).  It  can  be  checked  that 

<lQ  o  2  =  2  o  4Q. 

9  * 


(6.6) 


(6.7) 


Define  the  pre-momentum  mapping 

J( x,  rj)  =  J  o  2 (x,  t?)  =  J(x,  Pq(x)).  (6.8) 

A  straightforward  argument  shows  that  the  pre-momentum  mapping  J  :  Q  X  G  —  G  E  Ad 
equivariant.  Thus  the  level  set  J  Q^u)  is  invariant  under  the  action  of  the  isotiopy  subgroup  GQ  • 
Furthermore,  let  H  :  Q  x  Q  -*  1R  be  defined  as  II  =  II  o  2 .  From  (6.7),  the  function  H  is  invariant 
under  the  group  action  T .  The  functional  IQ  can  be  now  written  as  IQ  —  H  —  (I,  ( ) •  By 
the  invariance  properties  of  H  and  <7,  the  restriction  of  IQ  on  J  1(p), 


II, 


€| 


=  H  J- 


*(#*) 


-  (m,  0- 
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is  invariant  under  tlie  group  action  of  G M .  As  a  consequence,  the  geometric  picture  is  the  same 
as  in  Figure  6.1.  An  algorithm  analogous  to  Algorithm  6.1  can  then  be  applied  to  check  if  (xe,f) 
is  a  local  minimizer  of  restricted  to  Before  doing  so,  we  introduce  a  few  notations. 

The  riemannian  metric  restricted  to  the  subspace  A fx  provides  an  a; -dependent  bilinear  form  on 
the  Lie  algebra  Q .  This,  in  turn,  induces  a  pairing  (locked  inertia  tensor  associated  to  x  E  Q), 
I lock(x)  :  G  —*■  G*,  defined  through 

{£,  Ilock(x)Tj)  =  <  £q(x),  T]q(x)  »r,  (6.9) 

for  £,  r]  E  G  •  From  the  symmetry  property  of  the  riemannian  metric,  we  have 

I (ocfc(;E)p)  =  (Focfc(x)£,  ??)) 


namely,  I(ocfc(x)  is  symmetric.  Also,  we  assume  that,  at  x ,  the  locked  inertia  tensor  has  an  inverse, 
I/oc/.(x)-1  :  G*  —*  G-  On  the  other  hand,  the  gyroscopic  field  also  induces  for  each  x  E  Q  an 
element  Iy(x)  in  Q*  defined  by, 

(I y(x),  r/)  =  <  Y(x),  tjq(x)  >x,  V  rj  E  Q.  (6.10) 


We  refer  to  Iy(x)  as  the  (x -dependent)  gyro-momentum.  The  function  can  now  be  expressed 
as 


#«(*»  77)  =  I /ocfc(*)(»7  -  0) 

+  F(x)  -  W(x)0  -  (Iy(x),  0, 

=  |  1 iock{x){rj  -  0)  +  Vc(x) 


with  the  pre-momentum  mapping,  from  (6.3),  (6.8),  for  77  E  G , 

(j{x,  77),  C)  =  (J(®,  %(*)),  0 

=  <  »?<}(*),  Cq(x)  »*  +  <  ^(x),  Cq(x)  >*- 
=  (Iiocfc(x)77,  C)  +  (Iy(x),  C}> 


(6.11) 


or  we  may  write 


T(x,  77)  =  Iiocfc(x)77  +  Iy(x). 


(6.12) 


For  f.i  &  G* ,  the  associated  isotropy  subalgebra  £/Me  is  defined  in  (2.4).  With  the  inner  product 
induced  on  Q  by  the  locked  inertia  tensor  at  xe ,  we  define  the  orthogonal  complement  of  G^  to 
be 

Gi  =  {  C  e  G  :  (C,  Iiocfc(xe)7?)  =  0,  V  77  E  Gnc  }.  (6.13) 
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(6.14) 


Following  the  notations  used  in  [39],  we  define  the  maps  A  :  G  G* ,  and  ,4  :  Q  — >■  Q,  by 

^4(t?)  =  ad*  [j.e,  A(r))  =  I/ocfc^e)-1^??), 

respectively.  As  proved  in  [39],  we  have  the  following  lemma. 

LEMMA  6.2 

Provided  that  is  finite  dimensional  or  A  is  elliptic  with  respect  to  the  inner  product 
induced  by  I iock(xe),  we  have 

(i)  A  maps  Q  onto  G^r . 

(ii)  A  maps  Q  onto  G%c  C  G* ,  where  (7“e  =  {  f-i  €  G*  '■  {/J,  7?)  =  0,  V  rj  6  (?Me  },  is  the 
anniliilator  of  . 

I 

With  these  notations,  we  are  ready  to  apply  Algorithm  6.1  to  check  if  (xe,£)  is  a  local 
minimizer  of  restricted  to  j-1(pt)-  We  proceed  as  follows. 

Step)  0.  Fix  £  £  G • 

Step*  1.  It  is  straightforward  to  derive 

DH$(x,  r/)  (Sx,  Sti ) 

=  DV(.(x)8x  +  (Sr],  Iiock(x)(v  ~  0)  +  \(T1  ~  £>  (Dliock(x)Sx)(r]  -  £))• 

The  relative  equilibrium  is  given  by  the  conditions  DV^(xe)  =  0,  T]e  =  £,  which  match  with 
the  conditions  we  obtained  in  Algorithm  5.4. 

Step  2. 

For  the  relative  equilibrium  determined  by  the  pair  (xe,  £),  we  have 

lie  =  «/(ze,  0  =  Iiocfc("^e)^  +  (6.15) 

Now  we  find  the  space  Ker  DJ(xe,  £).  From  (6.12), 

d 

DJ(x,  rj)(6x,  6tj)  =  — 

=  (Dliock{x)Sx)ri  +  liock(x)Si ]  +  DlY(x)Sx. 

Here  again  x  +  e  Sx  denotes  the  integral  curve  corresponding  to  the  tangent  vector  Sx  at  x .  For 
(Sx,  Sr])  to  be  in  Ker  D.J(xe,  f),  we  must  have,  from  (6.16), 


J(x  +  eSx,  rj  +  c  Si]) 


■■  =  0 


(6.16) 
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*7  =  -hockixeT1  {(DIlock(xe)6x)Z  +  DIy(xe)8x ), 

=  I;ocfc(a;e)_1  ident^(xe)fe. 

where  the  map  ident y:  Q  xTQ  ^  Q*  is  defined  by,  for  (x,8x)  E  TQ , 

ident  re(x)<fo  =  -  \^Dllock(x)8x)^  +  DlY{x)8^j .  (6.17) 

This  map  specializes  to  the  map  ident %  defined  in  [39]  when  Y  =  0,  i.e.,  for  simple  mechanical 
systems  with  symmetry.  The  properties  of  this  map  play  an  important  role  in  our  subsequent 
development.  We  need  the  following  lemma. 

LEMMA  6.3 

For  x  6  Q  and  (,  is,  i)  £  Q ,  we  have  the  following  identities, 


(C>  {DIiock(x)r}Q(x))v)  =  ([C,  r]\,  Iiock(x)i/)  +  {[v,  77],  I i0Ck(x)(), 
(C,  DIy(x)rjQ(x))  =  (I y(x),  [(,  r?]). 


(6.18) 

(6.19) 


Proof 


The  proof  of  (6.18)  can  be  found  in  [39].  Here  we  only  verify  (6.19).  By  definition  (6.10), 

d  | 

<Iy(exp  erj-x),  Q 


(C,  DIy(x)t]q(x))  = 


de 

d_ 

de 


e=0 


<C  Y (exp  eri  ■  x),  ^(expep  •  x)  >e, 


exp  €  77*37? 


I  e  =0 

=  lvq  <  y,  Cq  >  (*), 

=  <LVqY(x),  Cq(x)  >r  +  <Y(x),  LvqCq(x)  >x, 

by  using  Lemma  2.2.  Also  we  have  the  identity  =  [C,  rj\Q .  With  Lemma  2.1,  it  follows 

that 


<C,  DIy{x)t]q{x ))  =  <  Y(x),  [C,  t]}q{x)  »*  =  (Iy(cc),  [C,  77]). 


I 

We  now  evaluate  the  map  ident  F  (xe)  restricted  to  the  space  AfXe . 

LEMMA  6.4 

For  tj  £  Q,  at  relative  equilibrium  ( xe ,  £), 

ident  F5(xe)7]Q(a;e)  =  ad*  fie  +  Iiock(xe)[v,  £]■  (6.20) 
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Proof 


From  the  definition  (6.17),  for  arbitrary  v  £  Q , 

(ident  Y^(xe)r]Q(xe),  v)  =  —(u,  (DIiock(xe)VQ(xe))0  -  ( DlY(xe)r)q(xe ),  v), 


From  Lemma  6.3,  this  could  be  further  written  as 

—  ([v,  T?],  I loek(,xe)0  ~  {[£>  rl)i  I lock{xe)v)  ~~  (ly^e)?  W  i  ^j) 

=  (I|ocfc(*e)f +  Iy(*e),  [V,  v\)  +  (W(*e)K,  V],  v), 
—  t')  T  (I(ocfc(-^e)[^»  vlt  ^)  • 

where  the  formula  for  /ie  in  (6.15)  has  been  used.  We  thus  established  (6.20). 


The  discussions  in  Step  2.  can  be  summarized  by  writing 

Ker  Dj(xe,0  =  {  (&c,  v)  G  T(r.i0(Q  X  (7)  : 

J?  =  I(ocfc(a-‘e)"1  ident^(xe)  &  }. 


I 


(6.21) 


Step  3. 

As  seen  in  (6.21),  the  component  of  Q  in  Ker  DJ(xe,£)  is  determined  from  the  variation  6x 
in  TXeQ  •  We  thus  only  need  to  decompose  the  kernel  space  with  respect  to  TXc{G tXc  ■  xe).  Since 

=  TXc(G^-xe)  =  {  VQ(xe)  6  TXcQ  :  rj  £  },  (6.22) 

we  find  the  orthogonal  complement  of  Afxe  with  respect  to  the  riemannian  metric  as, 

V  =  {&£  TXeQ  :  <  Sx,  r)Q{xe )  ~>Xe  -  0,  V  77  G  }  (6.23) 

Consequently,  the  space  S  can  be  written  as 

S  =  {  (Sx,  7 ])  G  V  x  g  :  rj  =  ident  }}(xe)  Sx  }.  (6.24) 


and  we  obtain  Ker DJ(xe,0  =  S  ©  T(Xe^{G„e  ■  {xe, 0) ,  where,  with  respect  to  the  action  $ 
defined  in  (6.6), 

T( xe,i){Gtit  ‘  {xeiO)  ~  {  (C<9(x‘e))  a(h£)  •  }' 

can  be  shown  to  be  a  subspace  of  Ker  DJ(xe,£). 

Step  4. 

Now  we  check  the  definiteness  of  the  second  variation  of  on  the  space  S.  The  block 
diagonalization  techniques  prove  to  be  useful  in  this  context.  First,  we  note  that,  under  conditions 
specified  in  Lemma  6.6  (see  below),  the  space  V  can  be  decomposed  as 

V  =  Vrig  ©  Vint  1  (6.25) 
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where 

Vrig  =  {  (Q(xe)  :  C  G  },  (6.26) 

V/jVT  =  {  €  V  :  (C,  ident^(xe)&:)  =  0,  V  C  G  }■  (6.27) 

It  is  this  decomposition  that  the  block  diagonalization  is  based  on.  On  the  other  hand,  by  definitions 
(6.4),  (6.22)  and  (6.26),  we  have  Mx  =  Af£e  ©  Vriq-  The  relationship  between  these  spaces  is 
geometrically  depicted  in  Figure  6.2. 


Figure  6.2.  Decomposition  of  TXQ 


Next  we  check  the  second  variation  of  IJ^  given  by, 

D2He(xe,£)-(8x i,  Tj\)  •  (6x2,  1 72) 

=  <^7i,  I iock{xe)r]2)  +  D2Vi(xe )  ■  Sxi  ■  Sx2, 

=  (ident^(xe)fei,  I(0cfc(^e)_1  ident^Xe)^)  +  D2V^(xe)  ■  8x\  ■  Sx2, 

for  (Sx  1,  rji),  (8x 2,  772)  6  KerDJ(a;e,0-  For  convenience,  a  bilinear  form  on  TXeQ  x  TXcQ  is 
defined  as, 

Sx2)  =  (ident^(xe)&i,  Iiocfc(a;e)_1identi5(a;e)&2) 

(6.28) 

+  if-V^Xe)  ■  Sx1  ■  SX2. 

Accordingly, 

D2Hs(xe,0-{fa  1,  m)‘(Sx2,  V2)  =  8x2). 

We  have  the  following  key  proposition. 
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PROPOSITION  6.5 


For  r]Q(xe)  G  Vrkj,  and  8x  G  VinTi  &z(VQ(xe)i  8x)  — 

Proof 

We  first  find  the  second  variation  of  Vi.  By  the  property  that  V  is  G -invariant  and 
Lemma  6.3,  we  have 

DV^(x)  ■  t)q(x)  -  BV(x)  ■  t)q(x)  -  {DIiock(x)  ■  Vq(x))0  ~  (DIy{x)  *  Vq{x),  £)■> 

=  -  ([£,  77],  hock{x)ti  +  M*)>- 

It  is  then  easy  to  see  that,  cf.  (6.17), 

D2V^(x)  ■  tiq(x)  ■  8x  =  -([£,  r)},  (DIiock(x)-8x)£  + DIy(x)-Sx). 

=  ([£,  rj\,  ident^(x)&c).  (6-29) 

Next  we  evaluate  the  bilinear  form  on  VRIG  X  VINT  •  Combining  (6.28),  (6.29)  and  using  Lemma  6.4, 
we  obtain 

B€(t lQ(xe),  8x)  =(ident^(xe)  r]Q(xe),  I/ocfc^e)-1  ident^(xe)  Sx)  +  ([£,  v],  ident^(ie)  Sx), 

=  (ad;ne  +  hockiXe^V^^iockixeT1  ident  y((xe)8x)  +  ([^,i?],ident \(xe)Sx), 

=  (od;/te,I/ocfc(*e)"1  ident  y€(*e)&), 

=  (A(r]),ident\(xe)8x),  (6.30) 

where  A  is  defined  in  (6.14).  From  Lemma  6.2,  A(v)  €  •  For  Sx  G  Vint  ,  by  the  definition  of 

Vint ,  cf.  (6.27),  the  desired  property  follows. 

I 

With  this  proposition,  the  second  variation  of  on  S  at  relative  equilibrium  is  diagonalized 
into  two  blocks.  Checking  the  definiteness  of  D2H$  on  8  is  thus  equivalent  to  checking  the 
definiteness  of  Bc  on  the  spaces  of  VRIG  X  VB/G  and  Vint  X  Vrvt  independently,  under  the 
assumption  that  (6.25)  holds.  These  techniques  often  simplify  the  computations  quite  significantly. 
In  particular,  the  form  of  B?  on  VRIG  X  VB/G  can  be  worked  out  explicitly.  From  (6.30), 

Bj(7 lQ(xe),  Vq(xe))  =  ( A (r?),  identy€(*e)  r,Q{xe)), 

=  (A(rj),  ad^ne  +  Iiock(x e)\jli  8}) i  (6.31) 

=  (ad^ne,  I(OCfe(xa)_1  ad*ne)  +  («d>e,  <«*„£>• 
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This  is  the  Arnold  block  analogous  to  the  one  in  simple  mechanical  systems  with  symmetry  [41]. 
The  gyro-momentum  is  buried  in  fie  and  can  affect  definiteness  of  this  block.  Definiteness  of  this 
block  ensures  the  decomposition  (6.25)  of  the  space  V,  which  is  proved  in  the  following  lemma. 

LEMMA  6.6 

Positive  definiteness  of  on  Vrig  x  Vrig  implies  that  V  =  Vrjg  ©  Vint  ■ 

Proof 

The  proof  is  analogous  to  the  one  in  [39].  We  only  consider  here  the  finite  dimensional  case. 
Letting  ( q(x6)  £  Vrig  H  Vint  ,  we  have  (  £  G^e,  and 

(u,  identy4(a:e)CQ(Ze)}  =  0,  V  v  £-Q^.  (6.32) 

We  choose,  in  (6.32),  v  —  M(C)  £  ■,  which  is  ensured  by  Lemma  6.2.  By  comparing  with 

(6.31),  we  get 

B€(C«(*e),  C <?(*«))  =  0. 

Since,  by  assumption,  is  positive  definite,  this  implies  £  =  0.  Namely,  Vrig  H  Vint  —  {  0  }• 
On  the  other  hand,  dim  Vrig  +  dim  Vint  =  dimV.  Thus  the  decomposition  (6.25)  holds. 

B 

With  this  Lemma,  we  do  not  need  to  verify  the  decomposition  (6.25)  explicitly.  It  is  guaranteed 
by  checking  the  definiteness  of  the  Arnold  block.  We  summarize  the  discussion  in  this  step  in  the 
following  theorem. 

THEOREM  6.7 

If  the  bilinear  form  B^  is  positive  definite  on  both  Vrig  X  Vrig  and  Vint  X  Vint  ,  then  the 
relative  equilibrium  (xe,  £g(xe))  £  TQ  is  relatively  stable  modulo 

B 

Now  we  have  completed  the  process  of  Algorithm  6.1  of  determining  the  relative  stability 
for  a  gyroscopic  system  with  symmetry.  The  block  diagonalization  of  the  second  variation  of  H £ 
is  achieved  on  the  constrained  subspace  S .  A  few  explanatory  remarks  follow.  First,  we  note  a 
necessary  condition  for  relative  equilibrium.  Namely,  at  relative  equilibrium  (®e,  £), 

ad ^  /.Le  =  0.  (6.33) 

This  result  holds  in  the  general  setting  of  hamiltonian  systems  with  symmetry,  see  Proposition  1.2 
of  [39]  for  the  proof.  Next  we  consider  the  amended  potential  introduced  for  simple  mechanical 
systems  with  symmetry.  From  (6.12),  we  may  construct  a  mapping  from  Q  X  G*  to  Q  X  O’  as 
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(x,  fi)  t-r  (x,  hockix)-1  [fi  -  Iy(x)))  . 

With  this  transformation,  the  functional  on  the  space  QxQ*  can  be  expressed  as,  from  (6.11), 

h^x,  n)  =  v^x)  -  (p,  t), 

where  the  function 

vnix)  =  v(x)  +  *)>  I(ocfc(a:)_1  (m  —  Iy(*))j  (6.34) 

is  called  the  amended  potential.  It  can  be  shown  that,  at  relative  equilibrium  (xe,£),  we  have 

nvjxe)  8x  =  DV^(xe)6x, 

D2Vfl(xe)  •  Sxi  ■  6x2  —  B^(fex,  Sx 2). 

Thus  the  stability  conditions  in  Theorem  6.7  are  equivalent  to  the  conditions  for  the  relative 
equilibrium  to  be  a  constrained  strict  local  minimizer  of  the  function  V M .  This  conclusion  is 
analogous  in  spirit  to  the  Lagrange-Dirichlet  theorem  [6].  We  phrase  it  as  a  theorem. 

THEOREM  6.8 

For  gyroscopic  systems  with  symmetry,  the  components  of  relative  equilibria  in  the  config¬ 
uration  space  are  the  critical  points  of  the  function  .  If  the  configuration  component  xe  of  a 
relative  equilibrium  is  a  constrained  strict  local  minimizer  of  the  function  V u  (i.e.  by  taking  out 
the  neutral  directions  tangent  to  •  xe),  then  the  associated  relative  equilibrium  is  relatively 
stable  modulo  G ^  . 

I 

REMARK  6.9 

In  most  practical  problems,  the  augmented  potential  W  is  easier  to  compute  than  the 
amended  potential  VM.  From  Theorem  6.7  and  (6.28),  it  is  clear  that  positive-definiteness  of 
the  second  variation  of  on  V  is  sufficient  for  stability.  Following  arguments  similar  to  the 
discussion  regarding  ,  we  get  an  analogous  statement  as  in  Theorem  6.8  with  replaced  by 
the  augmented  potential  .  However,  sufficient  conditions  obtained  from  Vtl  are  clearly  weaker 
than  the  conditions  from  Vj . 

I 

REMARK  6.10 

Theorem  4.4  and  Theorem  6.8  can  be  combined  to  assess  the  stability  ol  closed-loop  relative 
equilibria  obtained  by  the  application  of  a  gyroscopic  feedback  law,  cf.  4.12.  to  a  simple  mechanical 
system  with  symmetry  and  exterior  force. 

I 
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In  the  following,  we  consider  two  special  cases.  First,  it  is  easy  to  see  that  for  the  case  of 
Q  =  G ,  cf.  Figure  6.2,  Vint  =  {  0  }.  Consequently,  we  only  need  to  consider  the  Arnold  block 
for  stability.  Secondly,  for  the  case  of  G  —  50(3) ,  we  have 

V  =  £  X  rj,  ad*£  jl  =  )i  x  £.  (6.35) 

where  £,  fj  £  so( 3),  and  fi  £  so*( 3).  Thus,  condition  (6.33)  implies 

pe  X  £  =  0,  or  /re  =  A£,  (6.36) 

where  A  €  It  is  a  scalar.  It  follows  that  is  the  subspace  spanned  by  the  vector  £,  which,  in 
turn,  implies  that  =  Qg.  Recall  that  from  (5.13),  is  invariant  along  the  group  orbit, 

Gi  •  From  Remark  6.9,  we  conclude  that  for  this  case,  the  function  defined  in  (5.14)  is 
sufficient  for  determining  stability.  We  summarize  the  discussion  in  the  following  corollary. 

COROLLARY  6.11 

We  consider  a  gyroscopic  system  with  symmetry  (Q,K,Y,V,G). 

(i)  For  the  case  that  Q  =  G ,  positive  definiteness  of  at  relative  equilibrium,  defined  in  (6.28), 
on  Vriq  X  Vrig  implies  relative  stability  modulo  G M . 

(ii)  For  the  case  that  G  =  50(3),  a  strict  local  minimizer  of  the  function  V\,  defined  in  (5.14) 
on  the  space  Q/G $  induced  by  the  augmented  potential  gives  rise  to  a  stable  relative 
equilibrium. 

H 

These  observations  are  very  useful  in  applying  the  energy-momentum  method  to  specific  problems. 

After  discussing  the  block-diagonalization  of  the  bilinear  form,  or  the  second  variation  of  , 
we  consider  the  decomposition  of  the  symplectic  structure  on  T*Q  side  for  gyroscopic'  systems 
with  symmetry.  We  follow  closely  the  derivations  in  [39].  Recall  that  for  gyroscopic  systems  with 
symmetry,  the  symplectic  structure  on  T*Q  side  is  unaffected  by  the  presence  of  the  gyroscopic 
term.,  and  is  hence  the  canonical  one.  The  associated  momentum  mapping  is  given  by  J  :  T*Q  — >  Q* 
with  (J(x,p),  £)  =  (p,  £q(x)).  We  define  the  fundamental  mechanical  connection  a  :  TQ  — *■  Q 
as 

a(vx)  =  I Tock(.x)J(xi  K'{yx))-  (6-37) 

It  can  be  shown  that  a  is  a  connection  on  the  G-bundle  Q  ^  Q/G  and  g:(£q(:c))  =  Given 
jx  £  G* ,  a  1-form  :  Q  T*Q  is  induced  through  the  connection  a , 

(ZM(x),  vx)  =  {p,  a(vx)).  (6.38) 
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Regarding  Z M  as  a  map  from  Q  to  T*Q ,  we  can  find  the  corresponding  tangent  map  TZ ^  :  TQ  — *■ 
TT*Q .  Recall  that  the  space  V  C  is  decomposed  into  Vrig  and  Vint  as  defined  in  (6.26), 
(6.27),  respectively,  by  assuming  that  the  Arnold’s  block  is  definite.  We  define 

SRIG  =  {TXcZ^-(q (xe):  (  G  }, 

WiNT  =  {  TXZ^  ■  8x  :  8xe  VINT  },  (6-39) 

yViNT  =  {  (0.  fa)  ■  (fa  VqM)  =  0,  V  77  G  G}. 

It  can  be  verified  that  S  =  Srig  ©  Wjjvt  ©  V^int-  Alscb  we  have  tlie  following  proposition. 

PROPOSITION  6.12 

At  the  relative  equilibrium,  ze  =  (a;e,  pe )  G  T*Q ,  we  have 

(1)  for  C,  v  G  Gpt , 

Wo (ze)(TXcZfle  •  Cg(xe),  TXcZPe  •  r'Q(a'e))  —  — (/.ie,  [<,\  r'’]), 

(2)  for  (  e  ,  (&c,  <5p)  G  V^iNT  ©  W/jvr> 

W’o (Ze')(TXeZfle  -C (^b  ^P))  ~(Pei  [C>  Oi(&c)]), 

(3)  for  &x,  &2  G  Vint  > 

w'o (ze)(TXcZfJ,i  •  6x\,TXeZllt  ■  8xo)  —  dZPt(ftx2,  8x\)  —  (-ZMe,  [&2,  ^*i]), 

(4)  for  T]  £  Q ,  (0,  <5p)  G  Wjjyy,  wo(ze)(TxeZfle  -pQ(xe),  (0,  <5p))  =  0, 

(5)  for  (0,  ^Pi),  (0,  8p2)  G  wo(ze)((0,  8>Pi),  (0,  8p2))  —  0. 

a 

The  proof  is  omitted  here.  It  is  very  similar  to  that  in  proving  Proposition  IV.4,  pp.  61-62  in  [39]. 
Define  the  map 

p  '■  x  ~^INT  x  '  xe]A 

where  [Q  ■  xe]A  is  the  annihilator  of  Q  ■  xe  in  T*Q  ,  as 

p{ C,  &,  5p)  =  TXeZ^-(Q(xe)  +  TXeZ^-6x  +  (0,  fip). 

This  map  induces  a  bilinear  form 

LOp{((  1,  &El,  6pi),  (C'2,  8X0,  Sp2))  =  Wo(2e)(p(Cl)  ^xl)  fa)i  PiC’2,  8x2,  8p2)). 

From  Proposition  6.12,  this  bilinear  form  can  be  written  as 
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G^  Vint  [G  ■  xe]A 

r>±  /  Lie-Poisson  Bracket  •  Rigid-Internal  Coupling  •  0  \ 

y  u.,  I  ’  'l 


_ 

Vint 

- 

-  Rigid-Internal  Coupling  • 

Canonical  Symplectic  Structure 

[G-xe]A 

V 

0  •; 

plus  a  Magnetic  Term 

Figure  6.3.  Block-diagonalization  of  up 

wP((Ci,  &i,  Spi),  (C2,  fa 2,  <5p2)) 

=  “(Me,  [Cl,  C2])  -  (Me,  [Cl,  «(^2)])  +  (Me,  [C2,  a(&Cl)]) 

+  dZ^(5x 0,  Sx  1)  -  (.2^,  [&2,  Sxi])  +  {6p2,  8x1)  -  {Sp1,  6x2). 

This  shows  that  the  restriction  of  the  symplectic  structure  can  be  block-diagonalized  as  in 
Figure  6.3.  This  is  analogous  to  the  case  of  simple  mechanical  systems  with  symmetry  [39]. 
However,  the  gyroscopic  effects  enters  through  the  definition  of  pe ,  cf.  (6.15). 

Upon  the  completion  of  the  discussion  of  the  abstract  framework,  now  we  implement  the 
energy-momentum  method  in  more  detail  for  the  special  case  of  G  —  SO( 3).  This  physically 
corresponds  to  the  study  of  stability  properties  of  rotating  structures.  Through  the  isomorphism 
between  IR3  and  skew  symmetric  matrices  defined  in  (3.14),  we  define  the  locked  inertia  dyadic 

l°lo ck(x)  as 

(C,  hock(x)fi)  =  C-I lck(x)r),  (6.40) 


where  we  have  used  the  trace  pairing,  cf.  (3.16).  The  matrices  I/ocfc(a;)>  I°ock(x)  are  related  by  the 
following  formula.  For 


hi 

d 12 

hl\ 

1 

(h  2 

+  I33 

-IV2 

~h3  \ 

I12 

I 22 

I23 

I°ock(X)  ~  y 

' 

~  d\  2 

hi  +  h3 

-123  j 

\h  3 

123 

I33J 

V 

-Il3 

-1-23 

hi  + 122  J 

Also,  we  define 

IY(x)  =  1^), 

where  I y(x)  E  IR3 .  Namely, 

(Iy(*),  V)  =  I y{x)  •  '<)■ 

With  these  two  objects,  we  have  the  following  new  representations, 


(6.41) 


(6.42) 
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0  =  ILkMC  +  rY(x),  (6.43 a) 

V(  =  V(x)  -  -  W-5,  (6.436) 

ident^zjfo  =  ~(DI°ock(x)6x)  £  -  DIy(x)Sx ,  (6.43c) 


The  bilinear  form  defined  in  (6.28)  is  now 

Sxo)  =  ident^°(a:e)&i  •  I)’ocfc(a;e)_1ident^0(a;e)fe2 

+  D2  V^(xe)  ■  8x\  ■  6x2. 

The  Arnold  block  in  (6.31)  can  be  then  written  as,  cf.  (6.35), 

Ifocfc^e)  T  (o,d fjflei  ®^7j£) 

=  (He  X  7?)  ■  l°lock(xe)-1  (/ie  X  T?)  +  (fle  X  7?)  •  (T?  X  £) 

=  A2  ((XT?).  (i^W1  -  (£  x  77) 


(6.44) 


(6.45) 


It  is  thus  clear  that  for  the  Arnold  block,  we  need  to  check  the  definiteness  of  the  matrix 
I°oc k(xe)~1  ~  jl  along  all  directions  except  £.  Note  that  here  A  is  not  an  eigenvalue  of  the 
locked  inertia  dyadic  in  contrast  with  the  case  of  simple  mechanical  systems  with  symmetry.  The 
gyroscopic  field  afFects  A  through  the  gyro-momentum  term,  cf.  (6.36).  The  above  formulae  will 
be  used  in  the  following  section  for  the  example  of  two  rigid' bodies  with  rotors  coupled  via  a 
ball-in-socket  joint. 


7.  Two  Coupled  Rigid  Bodies  with  Internal  Rotors 

We  now  apply  the  energy-momentum  method  developed  in  Section  6  to  a  multibody  analog 
of  the  dual-spin  problem.  In  [48],  we  show  that,  with  appropriate  damping  mechanism,  the  system 
depicted  in  Figure  7.1  asymptotically  approaches  one  of  the  stable  relative  equilibria  corresponding 
to  an  associated  gyroscopic  system  with  symmetry.  Here  we  will  compute  certain  stable  relative 
equilibria. 

The  system  under  consideration  consists  of  two  rigid  bodies  connected  by  a  three- degree- 
of-freedom  spherical  (ball-in-socket)  joint  and  three  symmetric  rotors  mounted  on  the  center  of 
mass  of  one  body  along  its  three  principal  axes,  see  Figure  7.1.  These  rotors,  called  driven  rotors , 
are  set  in  constant  motion  relative  to  the  carrier  body.  We  assume  that  the  assembly  is  moving 
in  a  free  space.  For  simplicity,  the  inertial  reference  frame  is  placed  at  the  center  of  mass  of 
the  assembly.  This  corresponds  to  the  reduction  by  the  translational  invariance  of  the  system 
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as  discussed  in  [16],  [47].  Let  m1,  m2 ,  ms.,  i  =  1,2,3,  and  Ix ,  I2  ,  Is,  ,  i  =  1,2,3  denote 
tlie  masses  and  the  moments  of  inertia  of  body  1,  body  2  and  driven  rotors,  respectively.  Let 
£  =  m2{m\  +  mSl  +  ms 2  +  msJ  /  (mi  +  mSl  +  mS2  +  mS3  +  m2),  the  reduced  mass.  As  in 
the  classical  dual-spin  example  of  Section  4,  this  system  can  be  put  in  the  category  of  gyroscopic 
systems  with  symmetry  with  the  following  entities,  cf.  [46], 


Q  =  50(3)  x  50(3), 

K((Biu  1,  B2u2),  B2w2 ))  =  {uj  u%)  j1;)  , 

Y(BU  B2)  =  (Bm,  B2y2),  V(BU  B2)  =  0,  G  =  50(3), 

where, 
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3 

Ji  =  Ii  +  edjdi  +  y/  Is. , 

i—i 

3 

J2  =  I2  +  £  <$>d2  +  y  Ijr) . ,  J12  =  e dxBx  B2d2, 

i= 1 

and  the  components  yx ,  y2  of  the  gyroscopic  field  are  given  by  solving  the  equations 
Jl  2/1  +  J12  2/2  =  l  =  IS0,  Jf2  j/i  +  J2  2/2  =  0. 

The  associated  Lagrangian  is 

L(BX,  Qx,  B2,  fl2)  —  -  <  Oj,  Jifii  >  +  i  <  fj2,  J2f}2  > 

+  £  <  fli,  d\Bx  B2d2£l2  >  -f  <  /  >  . 

This  system  can  be  also  viewed  as  a  closed-loop  system  of  a  simple  mechanical  system  with 
gyroscopic  feedback  as  described  in  Example  4.5. 

Now  we  find  the  quantities  introduced  in  Section  6.  For  the  system  under  consideration,  the 
locked  inertia  dyadic,  cf.  (6.40),  can  be  determined  from 
<  £>  ^°iock(Bh  B2)rj  > 

=  <£q{BuB2),  fjQ(BuB2)> 

—  <  (BiJxBx  +  B2J2bT  +  BiJ^Bj  +  B2 J jo Bx)t)  >  . 

Thus  we  have 

I  fock(BuB2)  =  BXJXB^  +  B2J2Bj  +  BxJ12Bj  +  B2y2Bf. 

The  gyro-momentum  in  Q*  induced  by  the  gyroscopic  field  Y  can  be  shown,  cf.  definition  (6.42), 
to  be 

I °Y  =  Bx  l. 

Accordingly,  the  momentum  mapping  is 

V  =  I  iock(*)t  +  Iy(*) 

=  (^B\JXB^  T  B2J2B^  +  BiJx2bT  -f-  B2JX2B^[)Z  +  Bxl. 

The  augmented  potential  function  is 

V,{BuB2)  =  <£,  I?ocfc£>  -  <Bxl,  Z  >, 

—  ~ 2  ^  (BiJiBx  +  B2i2B2  +  BiJX2B2  +  B2jJ2Bx  )Z  > 

-  <  f,  -Si/  >, 

—  2  ^  (-/?iJi-Si  T  B232B2  T  6  f?x d\ B2 d2  T  f  B2d2B\dx)Z  Y 

-  <Z,  Bxl  >, 
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Now  we  apply  the  Principle  of  Symmetric  Criticality  to  find  the  conditions  for  re.ative 
equilibria.  The  first  variation  of  the  augmented  potential  is  derived  as  follows, 

DVz(Bi,  Bo)-{u iBlt  u2B2 )  =  <iBiJ!B[^  +  £Byl  +  e  B^d^B^do^,  ux  > 

+  <t=B2J2B^  +  elhdziBidiZ,  u2  >  . 

From  the  above  formula,  we  immediately  read  out  the  conditions  for  the  configuration  components 
of  the  relative  equilibrium  (5ie,J52e)  as  satisfying 

f  X  (BiJi Bit  +  Bil)  +  eBxdi  x  (f  x  ( B2d2  x  f))  =  0,  (7.2a) 

f  x  (B2J2B2O  +  eB2d2  x  (f  x  (B\d\  x  £))  =  0,  (7.2 6) 

These  are  very  similar  to  the  conditions  derived  in  [47],  except  that  a  gyroscopic  term  enters.  By 
taking  dot  product  with  f  on  both  side  of  (7.2a),  and  letting  Si  =  B\d\,  so  =  Bo  do  ,  we  obtain 
the  coplanarity  condition ,  cf.  [47], 

f  •  («!  X  So)  =  0.  (7.3) 

Accordingly,  the  gyroscopic  term  does  not  affect  the  coplanarity  condition  for  the  relative  equilib¬ 
rium  for  this  problem.  With  this  condition  (7.3),  equations  (7.2)  may  be  re-expressed  as 

£  x(J91J151Te  +  -Bi0  -  e{B1d1  -0(B2d2  xf)  =  0,  (7.4a) 

£  x  (52J2jB2tO  -  e  (B2do_  ■  OiB^i  x  0  =  0.  (7.46) 

Now  we  find  a  particular  relative  equilibrium  for  this  problem.  Let  {ei,e2,e3},  {fi,f2,f3}  be  the 
coordinate  frames  corresponding  to  body  1,  body  2,  respectively,  such  that 

Jie;  —  J2l'  ej ,  J2fi  —  ^2i  ff  >  f  —  1,2,3. 

It  can  be  checked  that  if  the  following  conditions  hold, 

£  =  |f|  B1  ex  =  |f |  Bo  flt  (7.5a) 

l  —  l\  ,  ( *  -56) 

di  =  ai  eo,  do  =  a2  f2,  (7.5c) 

conditions  (7.4)  are  satisfied.  Thus  the  conditions  (7.5)  are  associated  to  a  relative  equilibrium 
( B\e,Boe ).  From  (7.5a),  we  know  that 

Ble  ei  =  Boe  fx.  (7.6) 
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By  substituting  (7.5a),  (7.5c)  in  the  coplanarity  condition  (7.3),  we  get 

6i  •  (ei  X  Bei2)  —  0, 

where  Be  —  B^eB2e.  With  (7.6),  this  only  happens  when  Be  f2  =  ±  e2.  Thus,  we  have  two 
sets  of  relative  equilibria  expressed  in  terms  of  the  relative  shape  variable  Be , 

Be  fi  =  ei,  Be  f2  =  e2,  Be  f3  =  e3,  (7.7 a) 

Beii  =  ei,  Be  f2  =  — e2,  Be  f3  =  -e3.  (7.76) 

In  the  following,  we  will  study  the  stability  property  of  the  relative  equilibrium  corresponding  to 
(7.7b)  with  (7.5).  This  configuration  is  depicted  in  Figure  7.2. 


Figure  7.2.  Relative  Equilibrium  Configuration. 


The  energy-momentum  method  is  adopted  here  to  determine  the  stability.  We  first  need  to 

compute  the  second  variation  of  the  augmented  potential.  It  can  be  found  as  follows, 

D “ ( B\ ,  Bo')  ■  (fir B\ ,  uoBy )  •  (ui B\ ,  uoBo) 

d 


de 


DVz(ee Ul  Bi  ,ecu*B2)-  («i e€  UlBuu2ee  B2 ) 


6=0 


<  iiiBiJiBf £  —  BiJ\Bi  iii£  +  uiB\l,  rht  > 

T  <  uoBo ioB2i  —  Bo  JoBTuo^  uo£  > 

+  e  <  (uiBidi)£Bodo€,  «i  >  +  £  <  {uoBodo)^B\d\^,  u2  > 

+  2s  <  (uiBidi)^,  (uoBodo)£  >  . 


(7.8) 


Define 
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A  r>T  A  nT 

UX  =  BleUi ,  Uo  =  B^Uo- 

The  components  of  u1;  U2  will  be  denoted  by  (u.n,  ui2,  U13)  and  (u2i,  u22,  u23),  respectively. 
Also,  we  will  use  the  notations, 

Ji  =  diag{  Jn,  J12,  J13  },  J2  —  diag{  J2i,  J22>  J23  }• 

At  relative  equilibrium  (Ble,B2e)  such  that  (7.7b),  (7.5)  hold,  we  can  further  write  the  second 
variation  of  the  augmented  potential  as 

D2  Vj(i?ie,  02  e)  ■  (-BleUi,  i?2eU2)  *  (01eul ,  02eU2) 

=  £  |^|2«1  02(^11  —  u21)“  +  ((Jll  —  7l3)|£|  +  ^l)  ]£lu12 

(7 -9) 

+  ((Jll  —  Jl2  +  £  «2  )  l-f  I  +  h)  |£|u13  +  (J21  —  J23)|£|“u22 

+  ( J2I  —  J22  +  £  alft2)|Ci”u23- 

From  Remark  6.9,  we  check  the  positive  definiteness  of  the  second  variation  of  the  augmented 
potential  on  the  space  of  V.  For  the  relative  equilibrium  under  investigation,  we  have  the 
momentum  mapping,  cf.  (7.1), 

yUe  =  ((J11  +  J21  +  2£  uia2)|£|  +  h)  BXe  ei. 

Thus  the  Lie  algebra  corresponding  to  the  isotropy  group  is  =  Span{  Bieex  },  with  the 

orthogonal  complement  with  respect  to  the  locked  inertia  tensor, 

-  Span{  B\e  e2,  0iee3  }. 


The  space  V  is  given  by, 

V  =  {  (0ieui,02eU2)  :  <  (5ieUl,52eU2)?  VB2e)  >  =  0,  V  TJ  G  G ^  }, 

=  {  (0ieui,  02eu2)  :  (Jn  +  £  ui 02)1111  +  ( J21  +  £  aia2)u2i  =  0  }. 

The  second  variation  of  the  augmented  potential  restricted  to  V  is  now 

•(J5ieUi,  -0 2 e U2 )  •  (01eul>  02eu2) 


D2Vs(Ble,B2e) 


VxV 


=  hl±^l  +11 

V  J11  +  £  axa2  J 
+  ((Jll  -  Jl3)|f|  +  ^l)l^lu12  +  ((Jll  -  J12  + £0102)1^1  +  il)l$|Ui3 

+  (J21  -  J23)|£|“u22  T  ( Jll  —  J22  +  £  ala2)|sl'”u23- 
Consequently,  we  can  read  off  the  sufficient  conditions  for  stability  from  (7.10)  as, 


.10) 
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(Jn  —  Ji3)|£]  +  h  >  0, 

( Jn  —  J12  +  £  «i«2)|£|  +  h.  >  0, 
J2i  —  J23  >  0, 
J21  —  J22  "1“  £  Ct\d2  '>  0. 


The  above  discussions  are  summarized  in  the  following  theorem. 

THEOREM  7.1 

For  the  multibody  dual-spin  problem,  conditions  (7.5),  (7.7b)  give  rise  to  a  relative  equilibrium 
(f?ie ,  i?2e) .  Furthermore,  assuming  that 

h  ,  h 

■j^j"  >  d  1 3  -  J11,  and  -mj-  >  di2  —  dn  —  £aia2, 

the  relative  equilibrium  (i?ie,i?2e)  is  stable  if 

do  1  —  J23  >  0,  and  J2i  -  d22  -f-  £  d\(i2  >  0. 


1 


REMARK  1.2 

It  may  be  checked  that  the  positive  definiteness  conditions  for  the  Arnold  block  are 


(Jil  -  di3  +  J21  -  J23)|£|  +  h  >  0, 

(dll  d  1 2  T  d 0 1  d22  T  2c  Cll(22)|^|  T  Zl  >  0. 


These  conditions  ensure  the  decomposition  of  the  space  V,  cf.  Lemma  6.6.  It  is  easy  to  see  that 
these  conditions  are  implied  by  the  conditions  in  Theorem  7.1.  However,  this  is  not  sufficient  for 
stability.  There  are  additional  conditions  coming  from  the  other  block.  Thus,  for  such  a  coupled 
system,  we  could  never  regard  the  system  as  a  whole  rigid  body.  The  coupling  effects  should  be 
suitably  accommodated. 

I 

Now  we  consider  the  other  relative  equilibrium  coming  from  (7.7a).  The  second  variation  of 
the  augmented  potential  corresponding  to  the  case  that  the  relative  shape  is  identity,  or  the  two 
bodies  are  folded  can  be  found  from  (7.8)  to  be,  cf.  (7.9), 
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D^(Ble,  i?2e)  •  (-BieU1?  i?2eU2)  '  (BieUi,  i?2eU2) 

=  “£  |£l'al  02(^11  -  U21)2  +  ((Jll  —  j  13)1^1  +  ^l)|£|Ui2 

+  ((Jll  -  Jl2  -  £  ai«2  )|<f  I  +  ^l)|Clu13  +  (J21  -  J23)|C|2u22 

+  (J21  -  J22  ~  £  a2 )l^|"li23 * 

Even  restricted  to  the  space  V,  there  is  always  one  negative  term.  This  fact  suggests  that  this 
relative  equilibrium  may  be  unstable,  irrespective  of  the  rotor  speed.  Further  analysis  is  needed  to 
justify  this  statement. 


8.  Conclusions 

The  Lagrange-d’Alembert  principle  is  a  starting  point  for  modeling  natural  mechanical 
systems  subject  to  exogenous  forces.  If  the  forces  are  determined  through  feedback  laws,  then 
the  structure  of  the  closed-loop  system  can  be  used  to  assess  stability  properties  of  the  system.  In 
the  present  paper,  using  an  intrinsic  formulation  of  the  Lagrange-d’Alembert  principle,  we  have 
identified  a  class  of  feedback  laws  that  lead  to  gyroscopic  systems  with  symmetry.  The  (energy- 
momentum)  block-diagonalization  theorem  for  simple  mechanical  systems  with  symmetry  has  been 
extended  to  gyroscopic  systems  with  symmetry,  or  the  closed-loop  system.  Working  consistently 
on  the  tangent  bundle  side,  we  establish  the  splitting  that  block-diagonalizes  the  second  variation 
of  the  energy-momentum  function  at  a  relative  equilibrium.  The  splitting  depends  on  a  quantity 
that  we  refer  to  as  the  gyro-momentum  which  can  be  computed  in  terms  of  the  given  gyroscopic 
vector  field. 

The  gyro-momentum  also  enters  the  stability  criteria.  From  the  viewpoint  of  this  paper, 
the  gyro-momentum  is  the  key  control  parameter  and  thus  it  is  possible,  using  the  methods  of 
this  paper,  to  determine  whether  a  specific  gyroscopic  feedback  law  is  a  stabilizing  feedback  law. 
This  is  illustrated  in  the  example  of  Section  7  on  two  coupled  rigid  bodies  with  internal  rotors. 
This  example  is  a  natural  generalization  of  the  single  rigid  body  dual-spin  problem  studied  by 
P.S.  Krishnaprasad  [23],  Sanchez  de  Alvarez  [2],  and  more  recently  in  the  collaborative  work  with 
Bloch  and  Marsden  [7].  Other  more  complicated  examples,  including  dual-spin  satellites  in  central 
gravitational  fields  and  with  flexible  attachments,  appear  in  the  dissertation  of  Wang  [46]. 

In  future  work,  we  plan  to  investigate  bifurcations  of  relative  equilibria  with  respect  to  the 
gyro-momentum.  Examples  of  this  appear  in  the  work  of  Krishnaprasad  and  Berenstein  [25]. 
Control  strategies  based  on  bifurcation  of  relative  equilibria  may  be  effective  in  a  variety  of  problems. 
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We  hope  to  discuss  these  and  other  aspects  in  a  later  paper. 
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